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ABSTRACT 

The 1988 STLE Annual Meeting held a two- 
session symposium on cavitation In fluid films. 

This paper Is an edited review of the current 
research as discussed by the Invited participants. 
Phenomena and experimental observations Included 
gaseous cavitation, vapor cavitation, and gas 
entrainment. Cavitation In flooded, starved, and 
dynamically loaded journal bearings, as well as 
squeeze-film dampers was reviewed. Observations 
of cavitation damage in bearings and the possi- 
bility of cavitation between parallel plates with 
microasperities were discussed. The session on 
theoretical modeling summarized the transcavity 
fluid transport process, meniscus motion, and 
geometry, or form, of the film during rupture and 
reformation. Performance effects were related to 
heat-transfer models in the cavitated region and 
hysteresis Influences on rotor-dynamic coeffi- 
cients. Contributors presented a number of 
cavitation algorithms together with solution 
procedures using the finite-difference and finite- 
element methods. Although Newtonian fluids were 
assumed In most of the discussions, the effect of 
non-Newtonian fluids on cavitation was also 
discussed, 

INTRODUCTION 

This paper Is a survey of the state-of-the- 
art In research in cavitating fluid-film bearings 
and seals based on a symposium held at the STLE 
Annual Conference In Cleveland, Ohio, in May 
1988. The purpose of this article is to discuss 
the fundamental issues, experimental observations, 
theoretical modeling techniques, and further 
research needed In this field. It is hoped that 
this paper can serve as a reference publication 
for both the experienced researcher and those who 
are becoming Interested in this area. Neverthe- 
less, this article Is not a general review paper, 
as it is largely based on the viewpoints of the 
participants in the symposium. Each participant 
provided a short position paper covering his 
presentation from which the editors have drawn 
freely. 


This article is divided into three parts. In 
the first part, we present experimental results to 
provide the reader with the observed phenomena as 
well as some of the pertinent measurements the 
second part of the paper deals with theoretical 
modeling techniques. Part III consists of the 
editor's summary. 

In the text that follows, different types of 
cavitation are described: 

Gaseous Cavitation 

The lubricating oil often contains some dis- 
solved air and/or other gases. It is, therefore, 
likely that these gases could come out of the 
solution If the fluid pressure falls below the 
saturation pressure of the dissolved gases (nearly 
atmospheric pressure). This type of cavitation 
is called gaseous cavitation. In journal bearings 
under steady-state conditions, one would generally 
expect gaseous cavitation to take place in the 
divergent region where the mating surfaces tend to 
pull the lubricant apart, thus creating the condi- 
tion in which the gases could emerge from the 
solution and give rise to "atmospheric pressure". 
Side entrainment is another manifestation of 
gaseous cavitation when bubbles at the ambient 
pressure of the sides are drawn into the negative 
pressure regions of the bearing. 

Vapor Cavitation 

If the local pressure of the fluid coincides 
with that of the vapor pressure, then the boiling 
of the fluid by pressure reduction causes the 
formation of a so-called vapor cavitation. This 
type of cavitation normally occurs in dynamically 
loaded bearings used In common machinery. 

Understanding the mechanism of the formation 
and collapse of cavitation bubbles is of funda- 
mental importance since It is known to affect the 
performance of the bearing. Furthermore, it is 
the vapor cavitation that is known to be harmful 
as the bubbles violently collapse against the 
bushing surfaces, resulting in cavitation damage. 
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PART I 


FUNDAMENTAL AND EXPERIMENTAL OBSERVATION 
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NOMENCLATURE 
b Pad width, m 

C Bearing radial clearance, m 

D Bearing diameter, m 

d Depth of oi 1 , m 

e Eccentricity, m 

F Force, N 

f Frequency of load (Hz), 

G Duty parameter, [pUb/W] 

h Film thickness, m 

h a Amplitude of oscillation, m 

L Length of bearing In z-direction, m 
N Shaft rotational frequency, rev/min 

P Unit loading (W/LD) 

p s Oil supply pressure, m 3 /s 

Q Volumetric oil flow, m^/s 

Q Q/ (0.5 UCL) 

R Bearing or journal radius, m 

t T 1 me , s 

U Linear velocity, m/s 

u Velocity of the fluid in the circumfer- 
ential direction, m/s 

H Load on bearing, N 

x Coordinate in circumferential direction, m 

y Coordinate in radial or normal direction, m 
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EXPERIMENTAL STUDIES IN CAVITATION 


H. HESHMAT 1 

The experiments cited here deal primarily 
with the flow and formation of oil streamlets In 
cavltated regions. 

FLOODED BEARINGS 


Consider Fig. l f which depicts a hydro- 
dynamic journal bearing. As shown, the fluid 
film ends at ©2 after which It cavitates; 
the space at 9 > 82 is filled with oil stream- 
lets, released gasses, ambient air, and possibly 
oil vapor. In the absence of hydrodynamic posi- 
tive pressures In the region 82 - 8 [, the amount 
of flow ( 82 ) leaving O 2 remains constant at 
all 8 > 02 - If one observes the cavitation 
region through the bearing transparent shell, one 
sees fluid film streamlets in the manner shown In 
view B of Fig. 1. A sample of such experimen- 
tally observed streamlets seen through a transpa- 
rent bearing shell is shown in Fig. 2 (cf. (1)). 
This figure reveals that the oil adhering to the 
stationary surface is in narrow strips. The 
fluid film emerging from a cavltated region 
beyond the physical end of the pad 8 ^ is shown 
in Fig. 3. The pattern of ridges and troughs in- 
dicates that next to the shaft Is a sublayer of 
fluid that completely covers the rotating shaft. 
This is unlike the stationary surface. The fluid 
as it emerges from the bearing pad does not leave 
the shaft but adheres to it even when the film is 
not restrained by an upper bearing surface but is 
open to the atmosphere. The streamlets thus will 
adhere to the shaft even when they cross an oil 
groove. The streamlets emerge at the top of the 
picture from a cavltated zone and cross an axial 
groove before entering the downstream pad. As 
seen, neither the axial distribution nor the Indi- 
vidual width of the streamlets changes in cros- 
sing the groove, thus indicating adhesion of the 
entire flow (Q 2 ) to the moving surface. 

The above discussion suggests that the shape 
of the fluid film in the cavltated region is In 
the form depicted in Fig. 4. The fluid adheres 
In thin strips to the stationary surface, but 
forms a continuous sublayer on the moving surface. 

This pattern prevails as long as there is an 
upper pad and a thin spacing gap between shaft 
and bearing surface. Upon termination of the pad, 
the tips detach, and the entire flow <Q 2 > adheres 
to the shaft in the manner shown in Fig. 3. This 
"detached pattern" also prevails in the oil 
groove, whose dimensions are two to three orders 
of magnitude larger than the thickness of the film 
layer and is thus equivalent to an open space. 

The detached pattern reveals that there is a very 
weak bond between moving streamlets and 
stationary surface, permitting side entrainment 
of gas to take place in cavitated regions. 


’Mechanical Technology, Inc., Latham, NY. 
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Fig. 1. - Cavitation in journal bearings. 
TEMPERATURE EFFECTS 


The next topic pertains to the temperatures 
of the hot layer (Q 2 >. Numerous tests have shown 
that the temperatures of the fluid film in the 
cavltated region are more or less constant. This 
Is true despite seemingly further energy dissipa- 
tion, which takes place due to viscous shear over 
the extent of the cavitated region. Figure 5 
shows, for five sets of experimental operating 
conditions, that the temperatures of the cavi- 
tated region remained constant even though this 
region extended over nearly 80°. A possible 
explanation for this phenomenon is that the 
entrainment of ambient air separates the fluid 
from the bearing surface and causes less heat 
conduction to the thermocouples; i.e., thermo- 
couple tips are exposed to gas rather than hot 
oil. A partial explanation has also been given 
In terms of released gases and oil vaporization, 
which requires high enthalpies. However, even in 
cases where the subambient pressure is not low 
enough to cause oil vaporization, a constant 
temperature has been observed. To this end, the 
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L/D = 1; N = 1000 rpm; P = 60 psi 
(C/R) = 2 x 1CT 3 ; P s = 0.5 psi 

2. - Flow of oil streamlets over groove. 


View A- A oj 

Fig. 4. - Oil streamlets in cavitated regions of 
Journal bearings. 
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Fig. 5. - Constant temperatures in cavitated re 
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shape of the fluid film as postulated in Fig. 4 
offers an explanation of the role of gas entrain- 
ment and constancy of temperature in the cavitated 
region. The lower part of Fig. 4 shows the trans- 
verse velocity distribution In the oil streamlets. 
As seen, most of the oil moves with a velocity 
very near the journal speed (U); the principal 
strain thus occurs near the stationary surface. 
Since the power dissipation is related to the 
stress a = co(du/dy) at the moving surface and 
since, here (du/dy) * 0, the losses would be neg- 
ligible. Consequently, there would be no rise in 
temperature . 

To show that oil generally adheres to the 
moving surface, a series of tests was conducted 
with a tapered-land thrust bearing with reverse 
rotation. The apparatus and the experimental pro- 
cedures are described by Heshmat (2). The test 
bearing was instrumented with pressure taps and 
thermocouples to monitor the fluid film pressure 
as well as temperature. Rotation was opposite to 
the customary direction so that only the flat por- 
tion was loaded. Figure 6 gives some of the 
results obtained with this series of tests. The 
noteworthy points pertinent to the cavitation are 
the following: 

1. The pressure tap (located at the center 
of the flat portion of the pad) response attests 
to the presence of a hydrodynamic film, not 
boundary lubrication over the pad at the onset. 

2. In the diverging portion of the film, the 
temperatures actually tend to drop (Fig. 6). 

This further supports the view that, in the 
diverging film with cavitation, temperature rise 
is negligible because of the negligible power 
losses mentioned earlier. Furthermore, the end 
flow and streamlets in the cavitation region - 
induced by the diverging wedge configuration - 
basically adhere to the moving surface indepen- 
dent of surface geometry (journal or flat collar). 



Fig. 6. - Temperature profiles In a flat, diver- 
gent bearing. 


STARVED BEARINGS 

An experimental study documenting the perfor 
mance of journal bearings for conditions ranging 
from full hydrodynamic lubrication to states of 
extreme lubricant starvation is described by 
Heshmat and Plnkus (3). The effects of an Incom- 
plete fluid film on the extent of the film and 
the formation of cavitation will be reviewed here 
as functions of the amount of lubricant available 
to the bearing. 

Figure 7 presents the hydrodynamics of a 
starved bearing vis-a-vis a full film or flooded 
condition. For a fixed speed and supply oil 
temperature, the deviations of a starved bearing 
from one operating with a full fluid film are as 
fol lows: 

1. The film starts later and terminates 
earlier, l.e, 0] > 0]p and 02 < ^2F> producing, 
in essence, something similar to a partial 
bearing - though the upstream, and sometimes also 
the downstream, boundary conditions are different 

2. The decrease in the fluid film extent in 
starved bearings is shown in Fig. 8 and is pretty 
close to the theoretical values reported else- 
where. The start and end of the film was ascer- 
tained by means of pressure probes which located 
the points of the commencement and termination of 
the hydrodynamic pressures. 

3. The behavior of the fluid film in the 
starvation zone (9g - 0^ ) was observed via a 
transparent bearing shell. Its appearance is 


0 = 0 
eH 



Fig. 7. - Full and starved fluid films. 
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Fig. 8. - Extent of fluid film in starved 
bearings . 

sketched In Fig. 9. As seen, the full film under 
starved conditions started at more or less a 
constant angle, 9i = ©i p and did not vary along 
z. The starved zone between the start of the pad 
and ©if underwent cyclic changes in 
appearance. Most of the time, the zone 0g = 6 = 
B] was filled with an amorphous mixture of oil, 
foam, and gas bubbles (Fig. 9(a)). However, 
every 0.5 to 1.0 s, a pulse, originating from the 
cavltated zone downstream (9 > 62), transformed 
this region Into a pattern of oil streamlets 
(Fig. 9(b)) similar to the conventional 
downstream cavitation. The pulse frequency was 
checked against the journal's rotational 
frequency, but no relation nor connection was 
found between them. However, the pulsating 
streamlets and their frequency were related to 
the degree of starvation. The higher the degree 
of starvation, the longer the period of pulse 
frequency. 

4. The upstream portion of the starved 
bearing e B ^ 9 £ 9] filled with an oil and gas 
mixture is shown schematically in Fig. 9(a). 
Although gaseous bubbles were visible through the 
transparent bearing, they disappeared in the film 
and had no effect on the streamlets' patterns in 
the cavitated zone. 

Based on these observations, the model 
provides a reasonably accurate prediction as 
given in the literature (1,4-6). 


a. Basic Film Pattern 




b. Pattern During Periodic Pulses 

Fig. 9. - Fluid film pattern In starved bearings. 


DYNAMICALLY LOADED BEARINGS 

Film rupture in dynamically loaded bearings 
has been experimentally observed (7-24) to occur 
in two general patterns: (1) film rupture always 

exists through the load cycle, and (2) occurrence 
and disappearance of the film rupture are repeated 
every load cycle. In the case of a nutating 
squeeze-film damper, the ruptured film travels 
circumferentially ahead of the minimum film 
thickness (25). 

Recent studies by Nakahara et al . (18), 

Nakai et al . (14), and Natsumeda and Someya (26), 
all of whom used special pressure transducers 
(semiconductor strain gauge type) mounted in the 
test journal and bearing, report the existence of 
a sharp negative pressure except under a lightly 
loaded condition. Negative pressure (tensile) up 
to -1.2 MPa, as shown in Fig. 10, was reported by 
Natsumeda and Someya et al. (26). It is concluded 
that the heavier the dynamic load, the greater the 
negative pressure. Such a large value has been 
scarcely reported by others. Dyer and Reason 
(11) reported tensile stress with a maximum value 
of 0.64 MPa. However, their experimental results 
showed that the magnitude of tensile stress 
decreases as the eccentricity ratio is increased. 
Nakahara et al . (18) concluded from their experi- 
ments that cavity pattern changes were caused by 
the cavity growth due to the diffusion of dis- 
solved gas in oil under constant speed. Further- 
more, as reported recently by Walton et al . (25), 
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Fig. 10. - Circumferential pressure distribution 
of journal bearing under dynamic load (26). 


cavity patterns changed with the change of cavity 
locations. Both subcavity pressure and radial 
flow resistance across the circumferential groove 
changed di scont 1 nuously (traveling pressure wave 
up to five times that of supplied pressure). 

Depending on the apparatus and test operating 
conditions, five or more modes of transitional 
phenomena with cavity change were categorized by 
the investigators. Various boundary conditions 
have been used to predict dynamically loaded bear- 
ing performance. These are Sommerfeld conditions 
for light load, Reynolds conditions for heavy 
load, separation boundary condition for detecting 
tensile stress in the oil-film region downstream 
of the minimum film thickness position, and 
others in order to suit the fulfillment of the 
experimental data pattern. While the agreement 
between analysis and test is of the right order 
of magnitude (depending on the test and boundary 
conditions assumed) in predicting the extent of 
the film rupture zone, analysis does not ade- 
quately reflect the noted variations in the 
boundary conditions observed and does not account 
for the change in the shape of the film rupture 
zone with increasing speed. 

In a review of experimental results described 
by Walton et al. (25), pressure variations In the 


damper grooves exhibited significant flow rever- 
sals (circumferential feed and drain grooves in 
the sealed squeeze film damper). Also, the 
grooves, instead of being full of oil, appeared 
to experience film rupture, which contradicts the 
boundary condition assumption that the groove acts 
as a large sink. The flow fluctuations in the 
damper oil grooves, while not fully explainable, 
do affect damper performance and various stages of 
cavitation formation. 

It Is believed that side entrainment of gas 
takes place In most of those tribological test 
apparatus and sets up various cavitation config- 
urations. These discrepancies could be attri- 
buted, at least partially, to the differences in 
the interpretation of the various observations. 


CONCLUDING REMARKS 

Based on the review of the literature and 
results of an experimental Investigation, the 
following remarks are made: 

1. Side entrainment of the surrounding gas 
to the bearing gap at the downstream position is 
conceivable. 

2. The shape of the fluid film in the cavi- 
tated region is in the form given by Fig. 4. The 
fluid adheres by means of weak thin strips to the 
stationary surface, but forms a continuous sub- 
layer on the moving surface. This pattern pre- 
vails as long as there is an upper pad. Upon 
termination of the pad, the tips detach, and the 
entire exit flow adheres to the moving surface. 

In this region there would be no rise in 
temperature . 

3. Negative pressure (tensile stress) has 
been detected experimentally on the order of 

1 . 2 MPa (absol ute) . 

4. To fully understand and explain a sharp 
dent of negative pressure, one needs to consider 
the lubricant as a two-phase liquid consisting of 
oil and bubbles of gas and to take into account 
the surface-dilation effect of the surface ten- 
sion together with tensile strength of the film 
in the governing equation of motion. 

5. Under dynamic conditions, other types of 
pressure distributions exist and further exper- 
imentation and exploration are needed to quantify 
the nature of the formation of film rupture. 

6. Careful pressure measurements (from moving 
and stationary surfaces) of both the statically 
and the dynamically loaded journal bearing 
(although it is a difficult task) should be 
complemented by film temperature measurements in 
order to elucidate some of the unresolved aspects 
of cavitation. 


8 




CAVITATION IN JOURNAL BEARINGS AND SQUEEZE FILM DAMPERS: EXPERIMENTAL 


B. JACOBSON 2 

EXPERIMENTAL INVESTIGATIONS 

In nonstationary bearings, where the minimum 
film thickness moves during vibration, vapor cavi- 
tation can appear. This type of cavitation Is a 
breakage of the oil or oil-bearing Interface 
caused by tensile stresses in the lubricant. It 
can appear anywhere in the oil film — not just at 
the outlet where the film thickness is increasing 
in the direction of the sliding speed. 

Vapor cavitation has recently been studied 
by Jacobson and Hamrock (27,28) and Brewe et al . 
(29-31) experimentally, theoretically, and numeri- 
cally. The experiments show that the cavitation 
is influenced by the vibration amplitude, the 
vibration frequency, and the sliding velocity in 
the bearing. The clear Influence of the surface 
energy of the bearing surfaces Is also apparent; 
cavitation seems to start at the bearing surface, 
if it is made of PTFE (polytetrafluoroethylene) „ 
or in the middle of the film, if the surface is 
made of steel or PMMA (polymethylmethacrylate). 

In experiments vapor cavitation occurs when 
the tensile stress applied to the oil exceeds the 
tensile strength of the oil or the binding of the 
oil to the surface. The physical situation neces- 
sary for vapor cavitation to occur is a squeezing 
motion within a bearing. This Implies that the 
minimum film thickness is constantly changing in 
size and location as a result of the motion 
tangential and normal to the surfaces. 

Experiments showing vapor cavitation were 
first made at Lulea Technical University in the 
Machine Elements Laboratory and later at NASA 
Lewis Research Center, 

The geometry of the two test setups was the 
same, but the NASA Lewis test apparatus was 
stiffer and more precise, both in the measures of 
the bearings and also in the motion of the PMMA 
tube and the shafts. Figure 11 shows the NASA 
test apparatus with the high-speed camera and the 
driving motors for the rotation and vibration 
motions. The high-speed motion picture camera 
could take up to 11 000 frames per second with 
ordinary 16-mm motion picture film. In the 
experiments here, the framing rates were 1000 and 
2000 frames per second. 

The test bearing was lubricated by gravity 
feed from a can placed approximately 600 mm above 
the bearing. The test lubricant was a nondegassed 
Dexron II automatic transmission fluid with a vis- 
cosity of 0.066 N-s/m 2 under atmospheric condi- 
tions. The lubricant passed from the can, through 


2 Chalmers University of Technology, Goteborg, 
Sweden, and SKF Engineering and Research Centre, 
Nieuwegein, The Netherlands. 
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Fig. 11. - Test apparatus. 


a line and entered the bearing at the top left of 
the PMMA bearing (Fig. 11). On the right of the 
test bearing is the motor that gives the sliding 
velocity to the roller, and on the left is the 
motor that gives the squeeze velocity to the 
bearing. The speeds of these identical 0.5 kH 
motors can vary continuously from 0 to 250 r ad/s. 

Figure 12 shows the PMMA bearing, which can 
vibrate in a circular motion keeping the axis of 
vibration parallel to the axis of rotation. A 
special mechanism was manufactured to keep the 
centerlines of the bearing and shaft parallel and 
is also shown in this figure. 

The dynamic eccentricity mechanism can be 
adjusted to give dynamic eccentricity ratios 
between 0 and 1 or In actual measures from 0 to 
0.5 mm. The static eccentricity ratio is set by 
moving the motor with the shaft horizontally so 
that the centerlines of the two motors are offset 



Fig. 12. - PMMA- tube and eccentric. 
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Fig. 13. - Bearing geometries at four different 
times; e s * 0.4; e ( j * 0.6 


The static eccentricity can be between 0 and 0.5 
mm. The offset of the center of vibration from 
the center of the journal gives a static load, 
and the amplitude of vibration can be varied. 

Figure 13 shows four positions of the journal 
relative to the bearing for a dynamically loaded 
journal bearing. The static eccentricity ratio 
e s is set at 0.4, and the dynamic eccentricity 
ratio C(j Is set at 0 . 6 . 

EXPERIMENTAL RESULTS 

Figure 14 presents seven frames from the high- 
speed film showing the formation and collapse of 



(a) (b) (c) (d) (e) (0 (g) 

4 ms 12 ms 18 ms 20 ms 23 ms 24 ms 25 ms 


Fig. 14. - Dynamic cavitation bubble developed 
(C 5 ■ 0.34, C( j * 0.66, 035 - 19.5 r/s, and 0 ^ » 
92.7 r/s). 


vapor cavitation In a journal bearing. The experi- 
mental conditions In Fig. 14 are length to diame- 
ter ratio, 1/4; static eccentricity, 0.34; dynamic 
eccentricity, 0.66; rotational speed, to s , 19.5 
rad/s; and vibration speed, wtj, 92.7 rad/s. The 
journal surface, was made of PTFE, and the camera 
framing rate was 2000 frames per second. The time 
from the first frame to the last frame In Fig. 14 
Is less than 25 ms, Indicating how rapidly the 
vapor cavitation forms and collapses. 

Experiments with steel journals and with PTFE 
journals show a difference In the cavitation 
pattern. For steel journals cavitation bubbles 
start In different parts of the low pressure 
region, depending on where the largest bubble or 
dirt particle is located. For PTFE journals, the 
cavitation seems to start from the same position 
of the journal for each vibration, which indi- 
cates a rupture of the oi 1 -PTFE interface. 

The eccentricity ratios have a big influence 
on the size of the cavitation region and the dura- 
tion of this region from start to collapse. For 
small dynamic eccentricities (compared with the 
minimum oil film thickness), no vapor cavitation 
Is found. 
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CATEGORIZATION OF CAVITATION TYPES AND ONSET CONDITIONS 


D.W. PARKINS 3 

Oil films In dynamically loaded journal 
bearings are subject to both translatory and 
oscillatory squeeze motion of their containing 
surfaces. This discussion deals with oscillatory 
squeeze motion only. The principal aim of all 
work described here is to understand squeeze-film 
phenomena - simplified by removing the transla- 
tional component. It is important to distinguish 
sustained oscillatory squeeze motion from "single- 
shot movements" in which the oil film containing 
surfaces approach or separate. This paper deals 
with the former case. Some published work (32,33) 
deals with oscillatory motion in the frequency 
range of 1 to 2 Hz and only the transient first 
few cycles. Other works consider the radial 
squeeze motion of two circular surfaces, as, for 
example, in the application of squeeze-film dam- 
pers to the outer race of rolling-element bear- 
ings. This discussion deals with oscillatory 
normal motion of parallel surfaces at frequencies 
in the range likely to be met in machinery, say 5 
to 50 Hz, and for fully established behavior. 

Two of the important aspects of the more 
recently identified behavior are 

1. A steady force is generated by the oil 
film subject to normal oscillatory motion of its 
containing surfaces. The possibility was first 
suggested by Kauzlarich (34) who gave a theoret- 
ical treatment for a circular pad using a liquid 
lubricant and a harmonic displacement input. He 
concluded that inertial forces in the fluid cre- 
ated the steady force. Kauzlarich conducted a 
qualitative experiment to verify this finding. 
Subsequent experimental work by others (35-37) 
using nonharmonic oscillatory motion of a square 
plate has confirmed that such a steady force is 
generated, and - over a limited range of param- 
eters - determined the influence thereon of the 
initial oil film thickness, the oscillatory 
amplitude, and frequency of oscillation. Most 
recently, Kang (38) used a numerical solution to 
the Navier-Stokes equation to compute values of 
the steady force. 

2. Various regimes of cavitation appearing 
in the oil film due to the oscillatory squeeze 
motion have been identified and characterized, 
together with the location of the conditions that 
determine the onset of cavitation. Reference 39 
gives an initial classification and Ref. 36 fuller 
descriptions, all for behavior of an oil film con- 
tained in a square plate. Rodriques (40) predic- 
ted the growth and collapse of central cavitation 
patterns in a circular plate. His supporting 
experiments were conducted at only one frequency, 
approximately 1.0 Hz - but appear to show cavita- 
tion patterns which conform exactly to the charac- 
teristics of that which was later (39) defined as 
regime 2. 

3 Cranf i el d Institute of Technology, College 
of Manufacturing, Bedford, MK43 Oal , England. 


STEADY FORCE GENERATED BY THE OIL FILM 

Figure 15 contains data on the magnitude and 
direction of the steady force (F s ^) for eight 
cases and from three sources (35,36,41). Each 
source used the same basic equipment, but other- 
wise worked completely independently. Some curves 
show data obtained with invariant vibratory ampli- 
tude <h a ); one curve shows data obtained with 
invariant frequency; and one shows amplitude and 
frequency varying together. The data of Fig. 15 
are for Initial film thicknesses (y 0 ) from 0.11 
to 0.73 mm, vibratory amplitudes from 0.0055 to 
0.5 mm, frequencies from 5 to 50 Hz, three values 
of viscosity. Two values of oscillating surface 
area (A) are Included. The use of the steady 



Curve 

Ref. 

Amplitude 
of oscillation 
h a , 
mm 

Initial 

film 

thickness, 

Vo- 

mm 

Frequency, 

Hz 

Pod 

area, 

mm 2 

Lubricant 

kinematic 

viscosity, 

CS 

(a) 

35 

0.0055 

0.11 

10 to 50 

625 

313 

(b) 

35 

.0275 

.11 

10 to 50 

625 

313 

(c) 

41 

.044 

.22 

20 to 50 

2500 

750 

<d) 

41 

.044 to .66 

.22 

10 to 50 

2500 

313 

<•) 

36 

.4 

.73 

5 to 45 

2500 

370 

(1) 

36 

.4 

.5 

5 to 45 

2500 

370 

(g) 

36 

.4 

.25 

5 to 45 

2500 

370 

(h) 

36 

.07 to .47 

.25 

25 

2500 

370 


Fig. 15. - Measured steady oil film pressure. 



pressure P s t = F s t/A allows for its effect. 
Reference 42 gives additional values covering 
wider ranges. These include data for a third 
value of oscillatory surface area A and from a 
fourth source (43,44). It Is reassuring to 
observe that all the curves in Fig. 15 form a 
coherent pattern, even though they originate from 
separate sources. It is also reassuring to note 
that the additional data of Ref. 42 would also 
fit the coherent pattern if Fig. 15 were vastly 
extended. However, the data sparsely cover the 
parameter range. It is analogous to attempting 
to determine the picture of a jigsaw puzzle when 
only a few scattered parts are completed. 

Clearly, there remains the vital task of con- 
ducting a complete parametric study over a full 
range of variables . 

Data from Fig. 15 and Ref. 42 show that, 
depending on conditions, the steady force may 
either Increase or decrease the surface separa- 
tion. These data also show that the principal 
Influences on the steady force are initial film 
thickness (y 0 ), oscillatory amplitude (h a ), and 
frequency. The steady force varies directly 
with h a and 1 nversely wl th y 0 . Cavitation 
appeared in many of the cases reported. 

CAVITATION REGIMES 

It is widely established that cavitation may 
arise in such oscillating squeeze films. 

Rodriques (40) measured and computed cavitation 
patterns with a circular plate oscillating at 
1 Hz. May-Miller (39,41) identified three 
cavitating behavioral regimes, which he named 
regimes 1, 2, and 3. Each regime is 
characterized by a distinct method of feeding the 
bubbles and by the location of the bubbles in the 
displacement and pressure time histories. Exten- 
sion of the concept (36,45) shows that regime 3 
could take several forms. Figure 16 shows a 


FILM THICKNESS 




Fig. 16. - Typical oil film thickness and 
pressure-time history (with cavitation). 


typical displacement and pressure time history 
with cavitation of an oscillatory cycle for a 
harmonic force input. The pressure-time history 
Is not harmonic. It shows a rapid increase to 
maximum pressure followed by a fall to the mini- 
mum value. Regime 2 appears in this part of the 
cycle. A very rapid rise from the minimum pres- 
sure to near atmospheric is followed by a rela- 
tively long period near atmospheric pressure 
before the sharp rise again to maximum pressure. 
Regime 3 appears in this part of the cycle. The 
near atmospheric shoulder in the pressure curve 
of Fig. 16 shows a small positive gradient in the 
versions arising when regime 2 alone appears. 

The regimes and their subdivisions are 

Regime 1. No cavitation at any time in the 
oscl 1 latory cycle . 

Regime 2. Cavitation bubbles fed from within the 
film. Figure 17 is a photograph of a 
typical fully developed regime 2 
cavitation bubble pattern. The bubbles 
may contain dissolved gas extracted 
from the solution or vapor of the 
liquid, or a mixture of both. 

Regime 3. Cavitation bubbles fed from outside the 
film. Cavitation patterns in this 
regime may take several forms. 

(a) Edge fed: This form is characterized 

by smoothly curved bubbles entering 
simultaneously from all edges of the 
square plate. A typical example is 
shown in Fig. 18. Some bubbles therein 
show the smooth periphery of the bubble 
breaking into strlatlons which lead to 
the form described in under 3(b). 

(b) Leaf pattern: This form is character- 

ized by a leaf pattern supplied through 
a "trunk" that grows Inward from the 
plate edge (see Fig. 19). Usually one, 



Fig. 17. - Regime 2 cavitation pattern. 
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Fig. 18. - Regime 3(a) cavitation bubble pattern. 

but sometimes more such trees appear. 
Moreover, the tree may grow from a 
different side on different cycles of 
the same condition. References 36 and 
45 illustrate both of these points. 

(c) Leaf pattern with square parameter: The 

leaf pattern extends to cover all the 
oscillating plate surface except for a 
uniform border around three sides 

(see Fig. 20). 

(d) Stream fed: This form Is characterized 

by a leaf pattern whose perimeter 
follows that of the oscillating plate, 
but with a relatively wide (b/4) clear 
perimeter. This central leaf pattern Is 
fed by a number of diagonal streams. 
Figure 21 is a photograph of the fully 
developed pattern. 



Fig. 20. - Regime 3(c) cavitation bubble pattern. 



Fig. 21. - Regime 3(d) cavitation bubble pattern. 



Fig. 19. - Regime 3(b) cavitation bubble pattern. 


Reference 36 gives a sequence of photographs 
of cavitation patterns at a series of Identified 
points In the oscillating cycle for regimes 3(a), 
3(b), and 3(d). In each case the designated 
cavitation pattern exists for approximately 0.4 of 
the cycle time. The film is free of bubbles for 
a similar period. Bubbles advance and withdraw 
In the remaining time. With regime 3 (whatever 
subregime) each rapid increase in film pressure 
arises when the bubble front crosses the pressure 
transducer. The increase from minimum to 
atmospheric pressure occurs as the bubble front 
advances, and that from atmospheric to maximum 
pressure as the bubbles withdraw. An important 
finding in this work. (39,41) was the extremely 
rapid movement of the bubble front In both 
advance and withdrawal. 

Oil depth has a marked effect upon the onset 
cavitation patterns (36,45) . Figure 22, which 
Illustrates this effect, shows the value of 
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Frequency ( hz ) 


Fig. 22. - Effect of oil depth on onset of re- 
gime 3 cavitation <Y 0 = 0.7 mm). 

vibration amplitude (h a ) for a given frequency 
and oil depth (d) at which regime 3 occurs. Zones 
enclosed by dashed lines Indicate which subregime 
will appear. For example, if a frequency, say, 
of 20 Hz Is set and if amplitude is increased 
from zero, then, with d = 3 mm, regime 3(a) 
would appear at h a * 0.57 mm, regime 3(b) at 
h a = 0.75 mm, and regime 3(c) at h a * 1.3 mm. 

It is interesting to note that one bubble in 
Fig, 18 (regime 3(a)) shows its curved edge 
changing into embryo leaf patterns characteristic 
of regime 3(b). A similar exercise at a frequency 
of 10 Hz, d = 3 mm would show regime 3(a) cavita- 
tion at h a * 0.7 mm with no change therein 
occurring as amplitude h a is Increased 
further. Alternatively, for a frequency of 20 Hz 
and d =• 5 mm, Increasing amplitude h a from 
zero would show regime 3(b) cavitation at onset 
at h a « 1 .2 mm. 

Figure 16 shows how regimes 2 and 3 occupy 
different parts of the cycle. Consequently, it 


is possible for both regimes 2 and 3 to appear 
together in the same cycle. Indeed, Ref. 36 
gives an example of regimes 2 and 3(a) appearing 
together. Cavitation onset, of whichever regime, 
is Influenced by the geometry of the containing 
surfaces, A, y 0 , and vibration frequency and 
amplitude h a . However, evidence suggests that 
the onset of regime 3 is markedly affected by oil 
depth (d). The onset of regime 2 however, is 
influenced by the dissolved gas content in the 
liquid. These are Independent effects; therefore, 
it is possible to have many different combinations 
and sequences. It Is possible to transform from 
regime 1 (no cavitation) to any of the variants 
of regime 3 with, or without, regime 2 appearing. 
It is also possible to transform between the 
subdivisions of regime 3. Alternatively, it is 
possible to transform from regime 1 to regime 2 
without regime 3 appearing. 

This large number of possible forms of 
behavior make a confused picture when observed 
over a small range of variables, which renders it 
difficult to formulate classifications. There is 
a need for observation to be taken under a full 
range of all variables . 


CONCLUDING REMARKS 

It is now well established that an oil film 
subject to an oscillatory squeeze motion may gen- 
erate a force with steady and vibrating com- 
ponents. The steady force component has been 
measured separately and shown thereby to act - 
according to circumstances - either to increase 
or decrease the oil film thickness. Theoretical 
models that include fluid inertia have been shown 
to predict this steady force. Measurements have 
shown that the initial film thickness, vibratory 
amplitude, and frequency influence the magnitude 
of the steady force. 

Three cavitation regimes - including sub- 
regimes - have been identified and characterized 
according to the method of feeding the bubble. 

Each regime has a characteristic form of displace- 
ment and pressure-time history. Regimes 2 and 3 
may coexist In the same oscillatory cycle. A 
limited knowledge is available about the condi- 
tions which determine the onset of a particular 
regime or subregime and its subsequent trans- 
formation into another subregime. 

On both aspects, steady force and cavitation 
regimes, current knowledge is scanty, but suf- 
ficient to establish the existence of the listed 
features. More work needs to be done on both 
fronts before a full picture can be established. 
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CAVITATION IN MICROASPERITY LUBRICATION 


A.O. LEBECK 4 

It is well known that, even when sliders are 
made parallel and slide parallel in a lubricating 
film, fluid pressure load support develops so as 
to support some or all of the load and reduce 
friction. Lebeck (46,47) reviewed this behavior, 
the experimental data, and the mechanisms in 
detai 1 . 

A typical example of such behavior is shown 
in Fig, 23. The friction of carbon graphite pins 
sliding on tungsten carbide In water was measured 
using the apparatus shown In Fig. 24. The fric- 
tion coefficient at first decreases with Increas- 
ing duty parameter and then later increases. The 
decrease is caused by what are apparently hydro- 
dynamic pressures building up and reducing con- 
tact friction. The later Increase is, of course, 
caused by viscous friction. 

For many years it has been thought that 
microasperities serve as a source of fluid pres- 
sure load support in mechanical seals and parallel 
sliding lubricated contacts as reported in the 
experiments by Kojabashlan, Richardson (48), and 
Hamilton (49). Figure 25 shows how mi croasperi- 
t 1 es produce load support. Cavitation at the 
downstream side is essential to get a net load 
support. As Kojabashian and Richardson (48) show, 
one can predict that a theoretically significant 
load support will develop on a seal-like surface 
which has many microasperities. 

Experimental tests of this idea show that, 
if the ambient pressure is raised well above the 
cavitation pressure, no net load support results 
(Fig. 26). It follows that friction should 



^Mechanical Seal Technology Inc., Albuquer- 
que, NM. 



Fig. 24. - Hydrostatic bearing assembly diagram. 
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Fig. 26. - Zero load support. 
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Fig. 27. - Friction torque versus ambient pres- 
sure. Carbon pressure, 300 pst. Each point 
represents an average of 450 samples. 


increase because a lower fraction of the load is 
being supported by microasperities. This can be 
verified using the theory of Hamilton et al . (50). 
To test this hypothesis, Kanas (51) ran a series 
of experiments, using the apparatus shown in 
Fig. 24, where friction was carefully measured as 
a function of ambient pressure. His results (see 
Fig. 27) show no significant increase in friction 
with increasing ambient pressure. Similar results 
were obtained at higher speeds. 

Thus, In water, under the sliding conditions 
noted, and at a condition of (from Fig. 23) con- 
siderable fluid pressure load support, micro- 
asperity lubrication, or any other mechanism that 
relies on cavitation, does not appear to be a 
primary source of load support. 
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CAVITATION EROSION DAMAGE IN ENGINE BEARINGS 


F . A . MARTIN 5 

The experiences discussed in the 1979 paper, 
"Cavitation Erosion Damage in Engine Bearings: 
Theory and Practice, " by Garner, James, and 
Harriner (52) still reflect much of the current 
view on the subject. A few selected extracts 
from that p*aper are presented here. These 
Include examples of typical cavitation erosion 
damage, together with suggested mechanism of 
damage. Commonly applied palliatives and their 
effectiveness, based on engine experience, are 
also di scussed . 

''Cavitation erosion in engine bearings is a 
phenomenon which has assumed increasing signi- 
ficance during the past .. .years , probably as 
result of the design trends towards higher 
rotational speeds and. In some cases, higher 
rates of change of cylinder pressure... 

"...Damage is seen predominantly in diesel 
engine bearings, and on the rare occasions when 
damage occurs in gasoline engines it is normally 
as a result of operating under [incorrect 
conditions 3 . 

"In many cases the damage is restricted to 
local attack of the overlay and the performance 
of the bearing is virtually unaffected... 

However, in more extreme cases extensive loss of 
overlay and Interlayer material occurs, and the 
consequential adverse effects upon oil film cond- 
itions will reduce the service reliability and 
life of the bearings . . . 

"The fact that diesel engine bearings are 
more prone to cavitation erosion damage Is attri- 
buted to the more severe combustion conditions, 
and, additionally, to the necessity for more 
complex oil feeding arrangements with greater 
chance of flow discontinuities." 

Garner et al . state that cavitation erosion 
damage is associated with the collapse of vaporous 
cavities near to the bearing surface and suggest 
that the mechanism is the impingement by micro- 
jets of fluid. They describe a possible sequence 
of events in the collapse of vaporous cavities 
leading to the formation of high intensity micro- 
jets (see Fig. 28). Let us now look at the result 
of this in terms of actual damaged bearings. 

LARGE-END BEARING DAMAGE PREDOMINANTLY ASSOCIATED 
WITH OIL IN CONNECTING ROD DRILLING 

"The commonest form of damage occurring in 
large end bearings appears at first sight to be 
due more to flow erosion than to cavity collapse, 
but detailed examination of the damage confirms a 
cavitation effect. This can occur in either 
fully grooved or partially grooved bearings. 

5 Tri bol ogy Consultant to The Glacier Metal 
Company, Limited, Middlesex, England. 
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Fig. 28. - Mechanism of micro-jet formation (52). 

"In fully grooved bearings, the damage occurs 
predominantly at the groove edges, but can 
ultimately spread into the bearing lands. There 
is strong evidence that the damage is caused by 
inertia effects on the column of oil in the 
connecting rod drilling which is supplied from 
the large-end bearing groove to provide top end 
lubrication and piston cooling..." 

The damage may be cured by restricting the 
back flow from the rod drilling into the bearing 
(fitting a nonreturn valve has been one solution). 
In other cases, "the damage was diminished when 
the groove cross-sectional areas in both the main 
and large-end bearings were increased, thus 
allowing the oil column fluctuation to be better 
absorbed by the oil in the grooves.... 

"On partially grooved large end bearings 
damage frequently occurs at the groove run-out, 
either at the end of the partial groove, or at 
the side of the run-out, [Fig. 29]. The problem 
occurs at the end of the partial groove in the 
direction of shaft rotation, regardless of the 
position of the oil feed connection to the rod 
drilling. In the example shown in [Fig. 29], in 
which the damage had occurred within 400 hours, 


Fig. 29. - Cavitation damage associated with 
groove run-out In large-end bearing, together 
with effective solution (52). 
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the oil feed to the rod was adjacent to the 
damaged groove ending. The problem was overcome 
by the simple expedient of drilling a radial hole 
In the run-out section, as shown by the lower 
bearing In [Fig. 29] which had run for 1500 hours. 
Several other cases of similar damage have been 
successfully treated, ...either by provision of a 
radial hole or by terminating the groove 
conveniently at the end of a bridge-piece slot." 


MAIN BEARING DAMAGE PREDOMINANTLY ASSOCIATED WITH 
OIL IN CRANKSHAFT DRILLING 


"The majority of modern medium speed diesel 
engines use partially grooved main bearings in 
which the ungrooved arc is symmetrically placed 
in the bottom half, and Is approximately 120° in 
extent. Such bearings provide an acceptable oil 
film thickness within the narrow confines dictated 
by the overall design, but, at today's high ra- 
tings, are those most likely to suffer from cavi- 
tation erosion damage. This occurs downstream of 
but usually close to, the groove runout and 
typically takes the form shown in [Fig. 301. 

"In this Instance, the damage was found after only 
400 hours of operation, and was restricted to the 
soft overlay plated surface of the bearing, which 
was of the order of 0.03 mm thick. The stronger 
bearing lining underneath had not been attacked. 
This form of damage is characterized by a crescent 
form... In the particular example shown, after 
the initial short term damage, further progression 



Fig. 30. - Mild overlay damage in partially 
grooved main bearing (52). 



Fig. 31. - Cavity formation associated with Inter- 
action between journal drilling and partial 
groove (52). 


was slow. Bearings in this condition after 
several thousand hours operation may quite safely 
be refitted. However, In some Instances the 
damage can be much more severe... 

"It Is probable that as the oil drilling in 
the shaft traverses the groove run-out, cavities 
are formed at the intersection of groove run-out 
edge and the oil hole blend radius, as depicted 
diagrammatically in [Fig. 31]." 

"Various palliatives have been tried on 
different engine types to alleviate this type of 
damage. These have usually consisted of modifi- 
cations to the end of the partial groove. The 
most effective palliative so far adopted has been 
a tangential blend from groove to bearing surface 
with a square profile ([Fig. 32])." 



Fig. 32. - Design detail of successfully applied 
modification to groove ending (52). 
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Fig. 33. - Suction damage In large-end bearing as- 
sociated with large clearance (52). 


MAIN AND LARGE END BEARING DAMAGE PREDOMINANTLY 
ASSOCIATED WITH RAPID SHAFT MOTION 

Cavitation erosion damage of the type shown 
In [Fig. 33] (commonly referred to as suction 
cavitation erosion) occurs in the center of the 
lands In some fully grooved bearings. It is due 
to rapid shaft movement away from the bearing 
surface and is worse with large clearances. It 
occurs In the upper half of high-speed diesel 
main bearings (and also In a large end bearing 
during a series of development tests on a medium 
speed two stroke engine). 

Prediction of the journal orbit and exam- 
ination of the hydrodynamic conditions within the 
bearing helps to pinpoint the area of likely 
damage. It also enables the designer to quickly 
examine the relative merits and disadvantages of 
various designs. "For the purposes of balancing 
loads and integrated pressures, i.e., in pre- 
dicting the journal center orbit, it has been 
assumed, as Is common practice, that negative 
pressures cannot be sustained by the lubricant. 
However, It is considered that predicted negative 
pressures will indicate the position and Inten- 
sity of any cavity formation." 

This method was used to examine the hydro- 
dynamic conditions of a center main bearing 
showing the same type of damage as shown in 
Fig. 33. "The load diagram and resulting orbit 
are shown in [Fig. 34(a)], and It can be seen 


Polar load diagram Journal center locus 



Ruptured regions 

(a) 



Bearing angle 


Ruptured regions 
(b) 



0 180 360 


Bearing angle 

(c) 

Fig. 34. - Predicted hydrodynamic conditions asso- 
ciated with occurrence of cavity damage (52). 


that a rapid journal movement takes place between 
points B and C. Whilst experience suggests that 
rapid fluctuations In pressure are associated 
with such movements, the initial plot of rupture 
zones within the bearing ([Fig. 34(b)]) failed to 
show anything significant at the top of the 
bearing (0 or 360®) where damage occurred. It 
was only when the time step of the solution was 
refined (from 10® to 1° crank angle) that intense 
"negative" pressure zones were predicted 
([Fig. 34(c)])." The position of the maximum 
negative pressure corresponded to the position of 
cavitation damage on the bearing surface, thus 
Indicating that the cause of the damage was 
cavity formation (and subsequent collapse). The 
palliative here was to reduce the bearing clear- 
ance which alleviated the damage. However, it 
may be necessary in other instances to modify the 
polar load diagram. For example, the counter- 
weights could be changed. 
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CONCLUDING REMARKS 


Cavitation erosion damage In diesel engine 
bearings appears to fall broadly into three cate- 
gories, associated with 

1. Oil In crankshaft drillings and bearing 
oil feed features 

2. Effects of the oil within the connecting 
rod dr 1 1 lings 

3. Journal movements within the bearing 
clearance space leading to adverse oil film 
conditions. 


Relatively simple modifications to the 
detail design of the bearing and Its associated 
components can usually obviate problems from the 
first two categories, but often not from the 
third. However, these more Intractable problems 
do seem to show some correlation with predicted 
oil film conditions. A more detailed theory than 
that considered here will probably be needed 
before damage can be reliably predicted at the 
design stage. 


20 



REFERENCES 


1. Cole, J.A., and Hughes, C.J., "Oil Flow and 

Film Extent in Complete Journal Bearings," 
Proc . Inst . Mech . Eng . , 170, 499-510 (1956). 

2. Heshmat, H., "On the Mechanism of Operation 

of Flat Land 8earings," Interface Dynamics , 
D. Dowson, C.M. Taylor, M. Godet, and 
D. Berthe, eds., Elsevier, New York, 

149-160 (1988). 

3. Heshmat, H. f and Pinkus, 0., "Performance of 

Starved Journal Bearings with Oil Ring Lu- 
brication," J. Tribology , 107, 23-31 (1985). 

4. Articles, A., and Heshmat, H., "Analysis of 

Starved Journal Bearings Including Tempera- 
ture and Cavitation Effects," J. Tribology . 
107 , 1-13 (1985). 

5. Heshmat, H., "Theoretical and Experimental 

Investigation of Oil Ring Bearings," EPRI, 
Rotating Machinery Dynamics, Bearing and 
Seals Symposium, September (1986). 

6. Heshmat, H., and Pinkus, 0., "Mixing Inlet 

Temperatures In Hydrodynamic Bearings," 

J. Tribology , 108, 231-248 (1986). 

7. Jakobsson, 8., and Floberg, L. , "The Finite 

Journal Bearing, Considering Vaporization," 
Trans. Chalmers Unlv, Technology , Goteborg, 
Sweden, 190 (1957). 

8. Olsson, K., "Cavitation In Dynamically Loaded 

Bearings," Trans. Chalmers Unlv. Technology , 
Goteborg, Sweden, 308 (1965). 

9. Wakuri , Y., Tsuge, M., Nltta, Y., and 

Sanui , Y., "Oil Flow in Short-Bearing with 
Circumferential Grove," Bull. JSME , 16, 
441-446 (1973). 

10. Fall, C., "A Study of Cavitation, and the 

Phenomenon of Ribbing, in Steady, Hydro- 
dynamic Thln-Film Flow," Ph.D. Thesis, 

Leeds University, (1976). 

11. Dyer, D. and Reason, 8.R., "A Study of 

Tensile Stresses in a Journal-Bearing Oil 
Film," J. Mech. Eng. Sc 1 . , 46-52 ( 1 976). 

12. Savage, M.D., "Cavitation in Lubrication; 

Part 1. On Boundary Conditions, and Cavity- 
Fluid Interfaces; Part 2. Analysis of Wavy 
Interfaces," J . Fluid Mech . , 80, 743-768 
(1977). 

13. Pan, C.H.T., "An Improved Short Bearing Analy- 

sis for the Submerged Operation of Plain 
Journal Bearings and Squeeze-Film Dampers," 
J. Lubr. Technol . , 102, 320-332 (1980). 

14. Nakal , M., Kazamaki , T., and Hatake, T., 

"Study on Pressure Distribution In the 
Range of Negative Pressure in Sliding 


Bearings (1st report) -- The Relation 
Between the Behavior of Cavities and the 
Pressure Distribution Curve," J. Jpn. Soc. 
Lubr. Eng. , 27, 837-844 (1982). 

15. Pan, C.H.T., "Dynamic Analysis of Rupture in 

Thin Fluid Films. I - A Nonlnertial Theory," 
J. Lubr. Technol . , 105, 96-104 (1983). 

16. Parkins, D.W., and May-Miller R., "Cavitation 

In an Oscillatory Oil Squeeze Film," J. Tri - 
bology , 106 , 360-367 (1984). 

17. Kawase, T., and Someya, T., "An Investigation 

into the Oil Film Pressure Distribution in 
Dynamically Loaded Journal Bearing," 4th 
European Tribology Congress, Vol . II, 

5.2.3. (1985). 

18. Nakahar, T., Terasawa, T., and Aoki , H., 

"Transition PHenomena of Gaseous Cavitation 
in Thin Oil Film Betwee Non-Parallel Flat 
Surfaces," Proc. JSLE I~:ern. Tribology 
Conf . , Vol. 3, Elsevier, New York, 

1097-1102 (1985). 

19. Brewe, D.E., "Theoretical Modeling of the 

Vapor Cavitation In Dynamically Loaded 
Journal Bearings," NASA TM-87076/USAAVSCOM 
TR 85— C- 1 5 , and J. Tribology , 108, 628-637 
(1986). 

20. Hashlmoto, H., and Wada, S., "The Effects of 

Fluid Inertia Forces in Parallel Circular 
Squeeze Film Bearings Lubricated with 
Pseudo-Plastic Fluids," J. Tribology , 108, 
282-287 (1986). 

21. Elrod, H.G., "A Cavitation Algorithm, "J^ 

Lubr. Technol. , 103, 350-354 (1981). 

22. Orcutt, F.K., "An Investigation of the Oper- 

ation and Failure of Mechanical Face Seals," 
J. Lubr. Technol. , 91 , 713-725 (1969). 

23. Hughes, W.F., Winowlch, N.S., Birchak, M.J., 

and Kennedy, H.C., "Phase Change in Liquid 
Face Seals," J. Lubr. Technol . , 100, 74-80 
(1978). 

24. Findlay, J.A., "Cavitation in Mechanical Face 

Seals," J. Lubr. Technol . , 90, 356-364 
(1968). 

25. Walton, J.F., II, Walowit, J.A., Zorzi , E.S., 

and Schrand, J., "Experimental Observation 
of Cavitating Squeeze-Film Dampers," J^ Tri - 
bology , 109 , 290-295 (1987). 

26. Natsumeda, S., and Someya, T., "Negative Pres- 

sures In Statically and Dynamically Loaded 
Journal Bearings," Fluid Film Lubrication - 
Osborne Reynolds Centenary , D. Dowson, C.M. 
Taylor, M. Godet, and D. Berthe, eds., 
Elsevier, New York, 65-72 (1988). 


21 


27. Jacobson, B.O., and Hamrock, B.J., "High-Speed 

Motion Picture Camera Experiments of Cavita- 
tion In Dynamically Loaded Journal Bear- 
ings," J. Lubr. Technol. , 105, 446-452 and 
NASA TM-82798 (1983). 

28. Jacobson, B.O., and Hamrock, B.J., "Vapor 

Cavitation in Dynamically Loaded Journal 
Bearings," 2nd Intern. Conf. on Cavitation , 
Mechanical Engineering Publications, 

England, 113-140 and NASA TM-83366 (1983). 

29. Brewe, D.E., Hamrock, B.J., and Jacobson, 

B.O., "Theoretical and Experimental Com- 
parison of Vapor Cavitation in Dynamically 
Loaded Journal Bearings," Intern. Symp. on 
Cavitation, Sendai, Japan, Apr. 16-19 and 
NASA TM-87 1 2 1 /USAAVSCOM-TR-85-C- 1 9 (1986). 

30. Brewe, D.E., and Jacobson, B.O., "The Effect 

of Vibration Amplitude on Vapour Cavitation 
In Journal Bearings," Wear, 115, 63-73 and 
NASA TM-88826/USAAVSCOM-TR-86-C-26 (1987). 

31. Brewe, D.E., and Khonsarl, M.M., "Effect of 

Shaft Frequency on Cavitation in a Journal 
Bearing for Noncentered Circular Whirl," 
Tribology Trans. , 31, 54-60 (1980). 

32. Kuroda, S., and Horl, Y. , "An Experimental 

Study on Cavitation and Tensile Stress in a 
Squeeze Film," J. Jpn. Soc . Lubr. Eng. , 23, 
436-442 (1978). 

33. Hays, D.F., "Squeeze Films for Rectangular 

Plates," J. Basic Eng. , 85, 243-246 (1963). 

34. Kauzlarlch, J.J., "Hydraulic Squeeze 

Bearing," ASLE Trans. , 15, 37-44 (1972). 

35. Parkins, D.W., and Stanley, W.T., "Charac- 

teristics of an Oil Squeeze Film," J. Lubr. 
Technol . , 104, 497-503 (1982). 

36. Parkins, D.W., and Woollam, J.H., "Behavior 

of an Oscillating Oil Squeeze Film," J . 
Tribology , 108 , 639-644 (1986). 

37. Parkins, D.W., and Kahangwa, B., "Forces Gen- 

erated by an Oil Squeeze Film," ASME Paper 
80-C-2 Lub-39 , (1980). 

38. Kang, Byung-Hoo, "Prediction of an Oscilla- 

tory Oil Squeeze Film Behavior." M.Sc. 
Thesis, Cranfleld Institute of Technology, 
(1986) . 

39. Parkins, D.W., and May-Miller, R., "Cavitation 

in an Oscillatory 011 Squeeze Film," J^ Tri - 
bology , 106, 360-367 (1984). 

40. Rodriques, A.N., "An Analysis of Cavitation 

In a Circular Squeeze Film and Correlation 


with Experimental Results," Ph.D. Thesis, 
Cornell University, June (1970). 

41. May-Miller, R., "An Investigation of Cavita- 

tion in a Liquid Squeeze Film," MSc Thesis, 
Cranfleld Institute of Technology (1982). 

42. Parkins, D.W., "Oscl 1 1 atory Oi 1 Squeeze 

Films," Proc. Institute of Mechanical Engi- 
neers, London, Tribology - Friction, Lubri- 
cation and Wear - 50 Years On , Vol . 1. 
Mechanical Engineering Publications, 

Suffolk, England, 27-36 (1987). 

43. Stanley, W.T., "Dynamic Characteristics of an 

Oil Squeeze Film," MSc Thesis, Cranfleld 
Institute of Technology (1981). 

44. Kahangwa, B., "Squeeze-Film Forces in Oil 

Films," M.Sc. Thesis, Cranfleld Institute 
of Technology (1979) . 

45. Woollam, J.H., "An Analysis of the Vibratory 

Characteristics of an Oil Squeeze Film," 
M.Sc. Thesis, Cranfleld Institute of 
Technology (1984) . 

46. Lebeck, A.O., "Parallel Sliding Load Support 

in the Mixed Friction Regime, Part 1 - The 
Experimental Data," J. Tribology , 109, 
198-195 (1987). 

47. Lebeck, A.O., "Parallel Sliding Load Support 

in the Mixed Friction Regime, Part 2 - 
Evaluation of the Mechanisms," J. Tri - 
bology , 109, 196-206 (1987). 

48. Kojabashlan, C., and Richardson, H.H., "A 

Micropad Model for the Hydrodynamic Perfor- 
mance of Carbon Face Seals," Proc. Third 
Intern. Conf. on Fluid Sealing , Cambridge, 
England, BHRA, Paper E4 (1967). 

49. Hamilton, D.B., "Final Summary Report on the 

Dynamic Behavior of Liquid Lubricated Face 
Seals to the Rotary Shaft Seal Research 
Group," Battel le Memorial Institute, 
Columbus, OH, Feb. 22, 14 (1965). 

50. Hamilton, D.B., Walowit, J.A., and Allen, 

C.M., "A Theory of Lubrication by Micro- 
irregulari ties," J. Basic Eng, , 88, 

177-185 (1966). 

51. Kanas, P.W., "Microasperity Lubrication in a 

Boundary Lubricated Interface," Masters 
Thesis, Mechanical Engineering Dept., Univ. 
of New Mexico ( 1 984) . 

52. Garner, D.R., James, R.D., and Warriner, J.F., 

"Cavitation Erosion Damage in Engine Bear- 
ing: Theory and Practice," J. Eng. Power , 

102, 847-857 (1980). 


22 



PART II 


THEORETICAL MODELING AND PERFORMANCE 


23 




INTRODUCTION 

In part I, we were made aware of the numer- 
ous ways that cavitation could be manifested and 
of the categorization thereof. In particular, cav- 
itation can either be the result of (a) dissolved 
gas coming out of solution or (b) evaporation 
(flashing) of the fluid. Both types of cavitation 
are commonly observed In journal bearings, 
squeeze-film dampers, connecting rod bearings, and 
pistons. Both have a pronounced effect on their 
operation. The occurrence of cavitation in jour- 
nal bearings (1) results in reduced power loss, 
friction coefficient, bearing torque, and load 
capacity. Dowson and Taylor (2), in an excellent 
review of cavitation, point out that cavitation 
need not have a deleterious effect on the load- 
carrying capacity of bearings. Horsnell's calcula- 
tions (3) on load capacity for steady-state condi- 
tions show that, if cavitation is not considered, 
the load can Increase threefold over the predicted 
load . 

In recent years imposing higher loads and 
speed together with more complicated loading 
cycles resulted in bearings experiencing cavita- 
tion erosion damage, as discussed in the article 
by F.A. Martin in part I. In addition to the cavi- 
tation damage, there is the damage that can occur 
due to self-excited Instabilities that are encoun- 
tered under dynamic loading. These instabilities 
can be manifested as a whirling or whipping motion 
(4) of the journal center. Large vibrational 
amplitudes can result in large forces being trans- 
mitted to the system. A designer can circumvent 
that problem, or at least the severity of it, by 
knowing the speed threshold at which it occurs. 

This information is often obtained by referring to 
"stability maps" (5,6). Stability maps require 
the determination of hydrodynamic force terms cou- 
pled to the dynamical equations of motion. 

The calculation of the hydrodynamic force com- 
ponents is dependent on the film model used, espe- 
cially at high eccentricities. Many of the film 
models that are used to generate stability maps 
involve theories that are overly restrictive 
(e.g., narrow bearing theory, infinite bearing 
theory) and treat the cavitation in a very superfi- 
cial way. Notable of these is the ir-film cavita- 
tion theory in which positive pressures extend 
through half the circumference of the bearing and 
the other half Is regarded as cavitated. This 
boundary condition violates the conservation of 
mass. Specifying zero-pressure gradient and cavi- 
tation pressure at the boundary (i.e., Reynolds or 
Swift-Stieber boundary conditions) represents the 
conditions for film rupture in a reasonable way 
but does not properly represent the conditions 
when the film is reestablished. While this bound- 
ary condition may yield satisfactory results for a 
steadily loaded bearing, it leads to inaccuracies 
in situations where the bearing is dynamically 
loaded (7). In actual practice, the Swift-Stieber 
conditions are often compromised as a matter of 
convenience by neglecting the implementation of 
the zero-pressure gradient condition. This is 


much like the Gumbel conditions in which solutions 
were obtained to the Reynolds equation, and the 
cavitation was determined by disallowing the exist- 
ence of subambient pressures. However, dynamic 
loading causes changes In the local film thick- 
ness, which leads to nonstationary cavitation; 
that Is, the cavitation boundary Is in motion. 

This motion Is manifested as the growth and 
collapse of the bubble as well as downstream trans- 
port from the minimum film position. The appropri- 
ate boundary condition may require a condition 
other than a zero-pressure gradient at the bound- 
ary. The studies of Olsson (8,9) suggest that the 
usual Swift-Stieber film-rupture condition is ade- 
quate for dynamic situations If the cavitation 
boundary moves at a speed that is less than half 
the journal surface speed. It seems highly unlike- 
ly that this condition would be met during the Ini- 
tial stages of growth and the last stages of bubble 
col lapse. 

Jakobsson-Floberg (1) and Olsson (8) 
formulated boundary conditions for a moving bound- 
ary that conserved mass within the cavitated 
region as well as at the boundary (commonly 
referred to as the JFO cavitation theory). It was 
assumed that liquid was convected through the cavi- 
tated region in the form of striations extended to 
both surfaces In the film gap. There Is evidence 
that some liquid is transported in the form of an 
adhered layer to the faster surface (10,11). Pan 
(12) has broadened the JFO theory to accommodate 
those situations in which the mass transport 
through the cavitated region is not necessarily 
via liquid striations. For moderately to heavily 
loaded bearings (i.e., for a load-carrying 
capacity much greater than the surface tension 
forces of the lubricant), the adhered layer can be 
neglected (10). The striated flow in the JFO 
theory Is necessarily a Couette flow because of a 
constant-pressure assumption within the cavitated 
region. Jakobsson and Floberg's experimental find- 
ings (1) support the constant-pressure assumption. 

However, Etsion and Ludwig (11) have measured 
pressure variations of the order of 50 kPa inside 
a gaseous cavitation bubble. To resolve this 
apparent contradiction, one must know whether the 
cavitation bubble is gaseous or vaporous. Etsion 
(11) provides a plausible mechanism for determin- 
ing this, based on liberation and reabsorption 
rates between a gas and a liquid. The same argu- 
ments for pressure variation within a vapor bubble 
would not apply (13). In view of this discussion, 
the constant pressure assumption for vapor cavita- 
tion seems reasonable. Despite some of the contro- 
versy and certain lack of understanding, the JFO 
theory perhaps represents one of the best accounts 
of a dynamical theory to date of moderately to 
heavily loaded journal bearings and/or dampers. 

It is an improvement over the Swift-Stieber condi- 
tions, even for steady-state solutions, because it 
provides for film reformation as well as film 
rupture. Both rupture and reformation require a 
knowledge of the pressure gradient and fractional 
film content at the interface to determine its 
location. Unlike a rupture boundary, the reform- 
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atlon boundary Is subjected to a pressure flow 
( 1 . e . , nonzero pressure gradient). 

Furthermore, the fractional film content at 
the boundary Is determined by a residual fluid 
within the cavitated region that has been released 
at a rupture boundary (earlier in time) and gov- 
erned by the fluid transport law. Which condition 
(reformation or rupture) prevails depends on the 
relative motion of the boundary with the motion 
of the convected fluid normal to the boundary. 
Typically, under steady-state cavitating con- 
ditions, the rupture boundary occurs along the 
upstream boundary, and the reformation boundary 
occurs along the downstream end. Under dynamic 
loading and nonstationary cavitation, rupture 
or reformation conditions may be required at 
either the upstream or the downstream boundary. 

The programming task to effect these conditions 
in a computational algorithm is exceedingly 
tedious. In current practice this discourages 
their Implementation. 

The Elrod and Adams (14) computation scheme, 
which mimics the JFO theory, avoids the complex 
programming required to trace the moving boundary 


between grid points and to evaluate the required 
pressure derivative. Later, Elrod (15) modified 
it and presented It in more detail. Rowe and 
Chong (16) improvised the Elrod algorithm so that 
pressure is the dependent variable rather than 
fractional film content. Provision was made for 
pressure sources as well as automatic switching 
for full and cavitated regions. Miranda (17) was 
one of the first to apply this method for steadily 
loaded axial groove journal bearings. The full 
potential of the algorithm did not materialize 
until Brewe (7) applied it to nonstationary cavi- 
tation in a journal bearing undergoing noncentered 
circular whirl. For other bearing and seal appli- 
cations using the Elrod algorithm, the reader is 
referred to in (18-26) . 

This work represents a summary of the state- 
of-the-art and is by no means complete. Recently 
newer ideas and methodologies have emerged that 
contribute to the present methods and techniques 
of numerical modeling. We are sincerely grateful 
to the contributors of this composition in present- 
ing many of their innovations, ideas, and insights 
to this complicated problem. 
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NOMENCLATURE 


A] First R1 vl 1 n-Er 1 cksen tensor, 

A 2 Second R1 v 1 1 n-Er i cksen tensor, s -2 

C Bearing radial clearance, m 

Damping matrix, N*s/m 
D Bearing diameter, m 

e Eccentricity, m 

F Total fluid force, N 

F(s) Deformation gradient history 

f Dynamic fluid force, N 

f^ Adhered-fllm fill factor 

g Switch function 

H Nondlmensional film thickness, (h/C) 

H ^ Nondlmensional adhered film thickness 

h Film thickness, m 

he Film coeflcient of heat transfer for condens- 
ing oil vapor at low flow velocity, W/m- C 2 

I Last grid point 

1 Grid point Index 

Kij Stiffness matrix, N/m 

L Length of bearing In z-dlrectlon, m 

L' Axial halflength (L 1 - 1/2L), m 

l v Heat of vaporization, kJ/kg 

N Total number of time steps 

n Iteration level 

n Unit normal vector of the Internal boundary 

Ob Bearing center 

Oj Journal center 

P Nondlmensional pressure, p/viw(C/R ) 2 

p Fluid pressure, N/m 2 

Pvap Vapor pressure, N/m 2 

p a Atmospheric pressure, N/m 2 

Q Sweeping flux vector, m 2 /s 

Q 2 L Nondlmensional axial flow 


R Bearing or journal radius, m 
r Radius of flat bubble, m 
s Marking parameter In time, s 
T Temperature, °C 

AT v j Temperature change from oil vapor to 
liquid, °C 
T Stress tensor 

t Time, s 

U Surface velocity In x direction, and/or 

sliding velocity, m/s 
u Velocity of the fluid in the 

circumferential direction, m/s 

V Collapse velocity of the bubble, m/s 

v Velocity vector, m/s 

W Applied static load on bearing, N; width of 

cavitation bubble, m 

X Nondlmensional x-coordi nate , x/R 

x Coordinate in circumferential direction, m; 

Instantaneous position coordinate relative 
to journal center, m 

x s Vertical displacement of journal center rel- 
ative to bearing center, m 

x Instantaneous position coordinate relative 

to bearing center, m 

x Instantaneous x-component of journal veloc- 

ity, m/s 

Y Nondlmensional y-coordinate, y/R 

y Coordinate in radial or normal direction, m; 
Instantaneous position coordinate relative 
to journal center, m 

y Instantaneous y-component of journal veloc- 

ity, m/s 

y Instantaneous position coordinate relative 

to bearing center, m 

y s Horizontal displacement of journal center 

relative to bearing center, m 

Z Nondlmensional z-coordi nate , z/L' 

z Coordinate in axial direction, m 
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z c Location of bubble, m 

z e Axial location of the cavitation boundary, m 

z m 1/2 Bearing width, m 

a Inclination of unit normal (n) to direction 
of sliding, rad; non-Newtonian material 
parameter, N*s/m 2 

(3 Liquid bulk, modulus, N/m 2 ; non-Newtonian 

material parameter, N*s 2 /m 2 

Y Nondimensioned counterpart to the material 

parameter a, (aU/pR) 

5 Central difference operator; nondimensioned 

counterpart to the material parameter (3, 
( 6 BU 2 /uC 2 ) 

e Eccentricity ratio of bearing (e/C) 

£ Position of upstream cavitation boundary at 

any axial posi tion , m 

n Unwrapped circumferential coordinate along 

the cylindrical bearing surface, m 

0 Nondimenslonal density in full film, (p/p c ); 
fractional film content in cavitation region 

8 n Fill factor of adhered layer on sliding 
surface 

0B Bridged f 1 1 1 factor 

6q Fill factor of layer clung to stationary 
surface 

X Ratio of bearing half length to bearing ra- 
dius, <L’/Ro> 

p Viscosity, N*s/m 2 

p a Artificial viscosity, N*s/m 2 

5 Positional variable, m; Axial coordinate of 

the cylindrical bearing surface, m 

p Density of saturated oil vapor, kg/m 3 


x Period of dynamic load; the Instant the ad- 
hered film was formed at the breakup point; 
former instant in time, s 

x X y Shear stress, N/m 2 

<j> Attitude angle, rad. 

$ Filmflux, m 2 /s 

x Inverse function for the spatial position x, 
m 

xt Relative deformation function 

f Liquid phase transport flux, m 2 /s 

4/9 Circumferential transport flow, m 3 /s 

v|/ z Axial pressure flow, m 3 /s 

w Angular velocity of shaft, r ad/s 

Subscripts 
a Surface a 

b Surface b 

B Bearing 

c Cavi tation 

d Dynamic 

F Full fluid film 

j Journal 

max Maximum 

min Minimum 

S Starvation 

s Static, supply 

w Waist 
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DEVELOPMENT AND EVALUATION OF A CAVITATION ALGORITHM 


T.G. KEITH 1 AND D. VI JAYARAGHAVAN 1 

The effects of cavitation on performance have 
largely been disregarded In bearing performance 
calculations. For years, common practice was to 
modify full-film results In load calculations by 
setting all predicted negative pressures to the 
cavitation pressure (the half-Sommerfeld or GUmbel 
boundary conditions). Albeit this produced reason- 
able load values, It is also well known that the 
results violated the conservation of mass princi- 
ple. In more recent procedures the film pressure 
is assumed to develop from the point of maximum 
film thickness and end where the pressure deriva- 
tive vanishes (the Reynolds or Swift-Stieber bound- 
ary condition). Film rupture is more appropriately 
treated, but reformation of the film Is essentially 
neglected. 

Jacobsson and Floberg (1) and Olsson (7), pro- 
posed a set of self-consistent boundary conditions 
for cavitation to be applied to the Reynolds equa- 
tion. The boundary conditions properly account 
for mass conservation In the cavitation region. 

The collective procedure Is now generally termed 
JFO theory. This theory, when applied to moder- 
ately to heavily loaded journal bearings, produced 
values that were in satisfactory agreement with 
experimental values. The theory markedly Improved 
the prediction of cavitation regions in bearing 
analyses and for transient problems. 

Elrod and Adams (13) Introduced a computa- 
tional scheme that Incorporates JFO theory In a 
very simple manner. This procedure avoids the com- 
plexities of locating film rupture and reformation 
boundaries as It automatically predicts cavitation 
regions. Later Elrod (14) modified the finite dif- 
ference portion of the procedure to provide better 
results. The procedure, now generally termed the 
Elrod cavitation algorithm, has the following five 
objectives : 

1. To remove the pressure contribution to the 
mass flow In the cavitation region 

2. To observe the conservation of mass 

3. To set approximately correct boundary 
conditions at the boundaries of the cavitation 
region 

4. To produce a single equation for both the 
full film and the cavitated regions 

5. To Introduce the proper finite differenc- 
ing of the governing equation. 

The algorithm is based on a control volume 
formulation and uses the unit step (or switch) 
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function to eliminate the pressure term in the 
cavitation region. Elrod pointed out that the 
algorithm was developed "after considerable 
experimentation" and therefore does not present 
development details. Accordingly, it does not 
readily lend Itself to improvement. 

An attempt is made here to develop an algo- 
rithm that Is as effective as the Elrod algorithm 
In predicting cavitation in bearings. The pro- 
posed algorithm has the same objectives as the 
Elrod algorithm but Is developed without resorting 
to trial and error. Hopefully, this development 
will serve as a platform upon which further algo- 
rithm Improvement can occur. 

ANALYTICAL FORMULATION 


The two-dimensional, transient form of 
Reynolds equation for a Newtonian lubricant in lam 
inar flow, allowing for compressibility effects, 
can be written as 


/£hU . 

ph 3P\ 

a ( Eh: 3P\ 
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1 2y ax; 
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This Is an expression of the conservation of 
mass: the first term on the left is the squeeze 

term; the second Is the net mass flow rate In the 
x (flow) direction; and the third 1$ the net mass 
flow rate in the z (axial) direction. Note that 
the flow rate In the x direction consists of a con 
trlbutlon due to shear (Couette flow) and a contri 
butlon due to pressure gradient (Poiseuille flow); 
whereas, the flow rate in the z-directlon Is 
solely due to pressure. 


The density of the lubricant is related to 
the film pressure through the definition of the 
bulk modulus, 3, 1 .e. , 


3 = p(8p/a p ) [2] 

This relation suggests that it Is easier to work 
in terms of density in Eq. [1] than in terms of 
pressure. Accordingly, we can define a non- 
dimensional density variable 


8 = P - pc [3] 

where Pc Is the density of the lubricant at the 
cavitation pressure. 

For moderate to heavily loaded bearings, the 
pressure-induced flow within cavitated regions is 
small. Thus, the striated flow existing in these 
regions is primarily driven by shear, and the pres 
sure Is essentially constant. To reflect this fea 
ture, a so-called switch function is Introduced 
Into the pressure-density relation of Eq. [2], 

1 .e. , 

g3 - p(3p/3p) = 9(3p/3p) [4] 
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where 

I I in the full film region 

0 In the cavitated region 

Equation [4] can be directly Integrated to yield 

P = P c + gPOn 6) [5] 

Using Eqs. [1], [3], and [4] allows use to write a 
single expression for the conservation of mass in 
the bearing 



In order to reflect the correct physics of 
the flow, we desire a central difference of this 
term In the full-film portion of the bearing and a 
one-sided backwards difference in the cavitated 
region. This can be accomplished automatically by 
the Insertion of an artificial viscosity term Into 
Eq. [8]. This procedure is called type differenc- 
ing (27) when used for similar purposes in tran- 
sonic flow computations. Thus, Eq. [8] is written 
as 
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In the cavitated region Eqs. (5) and (6) are 
(p c hue) = 0 p =• p c 
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The variable 0 warrants some discussion, 
for, although it Is defined as the ratio of densi- 
ties, it carries a somewhat different interpreta- 
tion in the full film than in the cavitation 
region. In the former 0 > 1 and p > p c because 
of the compression of the lubricant. In cavitation 
regions 0 < 1 for another reason: The density of 

the lubricant is uniformly p c , but, because of 
the vaporization, the film is not completely occu- 
pied with lubricant. Thus, 0 is termed the frac- 
tional film content and 1 - 0 represents the 
void fraction. 


where |i a , the artificial viscosity switch func- 
tion, 1 s defined as 


1 0 in the full film region 

1 in the cavitated region 

The use of the switch function permits 
Eq. C 9 ] to be centrally differenced to yield 



Values of the variables on the half step may be 
computed from 


NUMERICAL FORMULATION 

Examination of the steady-state version of 
Eq. [6] for the full film (g = 1) reveals that it 
Is an elliptic partial-differential equation. 
Hence, the numerical scheme adopted for its solu- 
tion should reflect the fact that the dependent 
variable at any point in the full film depends on 
all neighboring variables. This permits the use 
of central differences. On the other hand, In the 
cavitated region (g => 0) , Eq. [61 is of the hyper- 
bolic type. Accordingly, downstream influences 
are not signalled upstream, and thus one-sided or 
upwind differencing should be used. 

In the following discussion, the shear flow 
and the pressure-induced flow components will be 
treated separately. A one-dimensional flow will 
be used in order not to obscure the details of the 
development. Extension to mul ti dimensional flow 
is not difficult. 

Shear Flow 

For simplicity the shear flow term will be 
written as 


( } U1/2 + ( [11] 

2 

Using Eq. [11], performing central differencing on 
the remaining derivatives in Eq. [10], and 
rearranging produce 



Sx 


P c h ®> 1 


aE u1 + b£- ♦ cE m 
2 Ax 


[ 12 ] 


where 


a = 1 - <n a >ui/2 
b = ( Ma>1+l/2 + Wi-1/2 
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c = -[1 + (jiaH-1/21 

Comparison of the definition of p a and the 
switch function introduced to cancel the pressure 
term in the cavitated region, via Eq. [4], reveals 
that 


9 - 1 - Ma U4] 

Thus, the coefficient in Eq. [13] are 
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a - < 91+1 + 9 1 > /2 


[ 21 ] 


b = 2 - (g i + i +2 g i + gi-p/2 V 


[15] 


c - -C2 - <gi + g i _ i ) / 2 ] 


Now, If all grid points are within the full film, 
all g 1 s are unity; therefore, a > 1, b > 0, and 
c = -1 . Hence, 
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which is the desired central difference form. On 
the other hand, if all grid points are within the 
cavitated region, all g's vanish so that a = 0, 
b = 2 and c - -2 . Hence, 


3_f h 3 ae\ 

This term exists In the full film region (6 > 1), 
but vanishes In cavitated regions <0 < 1). The 
switch function g adequately removes the term 
and, therefore, no artificial viscosity term need 
be added for this purpose. Also, when this term 
exists, the governing equation Is elliptic and 
thus can be centrally differenced. 

Equation [21] can be treated in a variety of 
ways which will ensure that it is central differ- 
enced In the full-film and made to vanish in the 
cavitated regions. However, depending on the man- 
ner In which the switch function is handled and 
the particular differencing scheme adopted, the 
results vary significantly at the cavitation bound- 
aries. Several alternative versions to the above 
were examined. The version that provided the most 
adequate and consistent results is one in which 
the derivative of g is not explicitly differ- 
enced. In particular, we may write that 
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and the desired upwind difference version Is 
obtained. 


The foregoing may be contrasted to Elrod's (15) 
shear flow expression, which is written as 
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If all grid points are within the full film, this 
yields 


_3e _3<e - 1) 3q(0 - 1) 

g 3F = 9 3x = 
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dX 


[ 22 ] 


The last term on the right hand side vanishes 
for all values of e because the spatial deriva- 
tive of g vanishes for all 9 except at 0=1; 
however, at this point, obviously, 0-1 is 
zero. Hence, we have 

80_ 3q (9 - 1) [231 

9 3x = 3x 


Now inserting Eq. [23] into Eq. [21] and central 
differencing yields 
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which, curiously, does not account for compressi- 
bility effects. In the cavitation region, 

Eq. [18] produces r 
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which agrees with Eq. [17]. 


Pressure Induced Flow 


The pressure Induced flow, written in terms 
of 0 Is 


which Is the same pressure term expression used by 
Elrod (15). 

SOLUTION PROCEDURE 


In this study cavitation in both slider and 
journal bearings was examined. All sliders were 
assumed to be one dimensional, whereas both finite 
and infinite length journal bearings were consid- 
ered. Because the procedures for solving one- and 
two- dimensional problems In these configurations 
differed, details of each will be briefly described 
in the following sections. 
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One-Dimensional Problems 


A purely Implicit finite-difference scheme 
was used to numerically determine the distribution 
of 9 for all one-dimens lonal problems. By this 
method, the finite difference version of Eq. [6], 
which Incorporates Eqs. [12], [15], and [24], can 
be written at each node in a general form as 


a A n+1 

a i e i-l 


b i 9 i 


. r A n+1 
+ c ri+i 


[25] 


where the coefficients a^ , b \ , q, and d^ are 
functions of Ax, At, 3, p, h, g, and the form of 
the boundary conditions applied. Consequently, 
the coefficients vary throughout the grid and 
change with time. Initially, all nodes are as- 
sumed to be In the full film, and at all nonboun- 
dary nodes 9 * 1 and g * 1. The switch function 
distribution was updated after every time step. 

Boundary conditions for slider bearings were 
taken to be of the Dlrichlet type In which 9 was 
set either greater than one (flooded) or less than 
one (starved) on the boundary nodes. Because the 
inlet and exit values of 9 were specified, the 
entire collection of nodal equations could be 
arranged Into the following tridiagonal matrix 
form: 
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be directly used in the solution. Brewe (7) has 
used the Gauss-Jordan procedure for this purpose, 
but found It to be extremely slow. To circumvent 
this problem, the two off-diagonal terms were 
lagged during the computations. By this procedure 
the following replacements were made: 


a Q n+1 
a i e i-i 


a i e i-i 


r « n+1 
C I-1 6 1 


C I-1 9 I 


C28] 


Thus, the terms were evaluated one time step be 
hind the actual computations and, therefore, are 
known during the <n+l) s * time step. Accordingly, 
they were Incorporated Into the right hand side of 
Eq. [27], and the Thomas algorithm applied. 


The time march process was continued until 
the change In 9j between two time Intervals was 
less than 10 ~ 4 . At this point, it was assumed 
that a converged steady state had been reached. 


Two-Dimensional Problems 


An alternating direction implicit (ADI) 
finite difference scheme (7) was used to determine 
the distribution of 9 for all finite journal 
bearing problems. By this method, the time step 
Is split Into two parts: in the first half, 

changes In the flow or x direction are calculated 
while old 0 values are used to determine the gra- 
dients In the axial or z direction; in the second 
half, the recently calculated 0 distribution is 
used to evaluate the flow direction gradients 
while the axial changes are updated. 


[26] 

In this form the collection of equations may be 
readily solved by application of the Thomas 
algorithm (27) . 

For journal bearings the boundary conditions 
depend on the means of lubricant supply. At oil 
grooves, values of 9 are assumed which corre- 
spond to the given supply pressure and are calcula- 
ted from Eq. [5], The procedure remains the same 
as for slider bearings. If, however, the bearing 
is submerged, periodic or continuous boundary con- 
ditions must be employed. Application of these 
conditions produce the following matrix form of 
nodal equations 
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The coefficient matrix In Eq. [27] is not tri- 
diagonal and therefore the Thomas algorithm cannot 


The Thomas algorithm can be used in each part 
of the time step. If again, the appropriate varia- 
bles are lagged for the submerged journal bearings 
case. In accordance with Brewe 1 s (7) suggestion, 
the switch function was updated after every half 
time step. Calculations were performed over half 
the length of bearing due to symmetry with resepct 
to the axial center. 

On occasion an oscillation of the 9 distri- 
bution was encountered. This oscillation was trig- 
gered at the cavitation node and has been reported 
elsewhere (e.g., (7)). To remedy this problem, a 
variable time step was Implemented according to 
the following scheme: 


At k - At m1n 


A W (k - 1)/<N - 1) 


At, 


min 


k = 1 , 2, . . . , N 
[29] 


where At m ^ n and At max are grid dependent and 
N Is normally taken to be 8. 

Not only did the variable time step help 
reduce the oscillation, It also reduced by approxi- 
mately 15 percent the number of time steps re- 
quired for the solution to reach steady state. 
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p - p c + (3(0 - 1 ) 


[30] 


Reference 28 examines the details of the cur- 
rent approach for both slider and journal bearings 
in one and two dimensions. The principal purpose 
of that effort was to compare the computed results 
with the published values of Elrod (14,15): The 

two algotlthms produced comparable results. A sec- 
ond purpose of that work was to examine a few modi- 
fications to the basic flow terms used In the Elrod 
algorithm. In particular, Elrod neglects In the 
shear flow term the compressibility effects In the 
full film region. Also, an approximate equation 
relating the pressure to the fractional film con- 
tent Is sometimes used. This approximate relation 
uses the first two terms of a series expansion of 
the logarithmic term In Eq. [5], l.e., 


CONCLUSIONS 

A cavitation algorithm that automatically 
alters the form of finite differencing used for 
the shear induced flow has been developed. The 
algorithm is an outgrowth of the Elrod cavitation 
algorithm. The two algorithms produce comparable 
results (28). However, the newer version may lend 
Itself more readily to Improvement. The effects 
of film compressibility, equation of state of the 
lubricant, and grid spacing near cavitation bounda- 
ries on the pressure distribution were Investi- 
gated (28) . 
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FUNDAMENTAL ISSUES CONCERNING THE MOVING BOUNDARY 


CODA H.T. PAN 2 AND JEAN FRENE 3 

This discussion deals with the motion of an 
Internal cavitation boundary of an Incomplete 
film, with emphasis on applications in which the 
dynamic load component overshadows the magnitude 
of the static load. In such applications one is 
concerned with supplying lubricant to the region 
where the surfaces are approaching each other, 
while voids should dominate the region of separat- 
ing surfaces such that viscous friction may be 
kept minimal. The commonly accepted model of the 
ir film is an attempt to approximate the extent of 
filling of the bearing gap. As engineers desire 
to refine the design art, a more precise descrip- 
tion of the Internal boundary of an Incomplete 
film in a bearing gap is sought. 

THE MODELING VIEWPOINT 

A predominant concern In the analysis of an 
Incomplete film is related to how flow continuity 
is maintained across the entire perimeter of the 
f i 1m, which compri ses : 

1. Inlet film edge or edges. 

2. Boundaries fixed by the geometry of the 
bearing. 

3. Internal boundary or boundries. 

Suppose, that at any Instant, such delinea- 
tion of the film perimeter can be asserted. Then 
one can assign the value of the film pressure at 
the film inlet to be that of the lubricant supply 
pressure; the geometrical boundaries can be repre- 
sented by the ambient pressure; and the "cavity" 
pressure can be assumed to prevail near the inter- 
nal boundary or boundries. Consistent with this 
set of pressure boundary conditions, a Dirichlet 
type problem Is prescribed by the lubrication 
theory and may be completely solved to yield, not 
only the film pressure distribution, but also the 
normal film flux at all points of the film 
perimeter. 

While the film continuity condition Is inher- 
ent In the interior of the bearing film as dic- 
tated by the lubrication theory, the transboundary 
continuity condition is more subtle. One may bor- 
row from gas dynamics the notion of the Rankine- 
Hugoniot conditions to allow for abrupt changes of 
state properties while preserving continuity of 
lubricant transport. As Rankl ne-Hugonlot jumps 
are approximations of the diffused transport proc- 
esses across a moving shock wave, the analogy for 
a moving cavitation boundary Is the change of flow 
structure that takes place within a distance of 
the scale of the bearing gap. 
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SOME EXPERIMENTAL FACTS 

Before attempting to examine the equivalent 
Rankl ne-Hugonlot problem at the Internal boundary 
of an incomplete film, It Is instructive to recall 
some representative experimental facts of incom- 
plete films. Two photographs (Fig. 1) which first 
appeared in the first leeds-Lyon symposium (29), 
reveal that the notion of a smooth internal bound- 
ary which separates a filled gap from a void re- 
gion is at best an approximation. In reality, the 
void region has a textured two-phase flow struc- 
ture. The border between the filled gap and a 
void region exhibits Irregularities associated 
with the two-phase flow structure. To circumvent 
such complicated details, it is convenient to deal 
with spatially averaged quantities of the two-phase 
region 1 n model 1 ng. 

It is clearly revealed that a portion of the 
void-dominated region is locally filled with liquid 
bridging both surfaces. Less visually evident, but 
not to be overlooked, are the possibilities that a 
thin liquid layer may be attached to either but 
not both surfaces. Thus, a complete morphological 
description of the void dominated region must in- 
clude the following distinct liquid states: 

State A - A thin layer of liquid that adheres 
to the sliding (journal) surface. 

State B - A liquid streak or spot that 
bridges both surfaces. 

State C - A thin layer of liquid that clings 
to the stationary surface. 

Excluded from present consideration are minute liq- 
uid droplets that move freely between the two sur- 
faces. 


O^GINAL. PAGE 
Aim J WHITE PHOTOGRAPH 



Fig. 1. - Appearance of high void content flow 
structures (29). (a) Divergent gap between 

rotating cylinder near stationary plane, (b) Sep- 
aration of sphere from submerged stationary plane. 
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EQUIVALENT RANKINE-HUGQNIOT 

The equivalent Rankl ne-Hugonlot Is concerned 
with the abrupt changes across a moving Internal 
boundary. It deals with the conservation of the 
liquid phase as related to the thln-fllm flow on 
one side, liquid-phase transport on the other 
side, and the motion of the boundary. 

Figure 2 shows that the normal unit vector of 
the Internal boundary at E Is directed Into the 
void region and is Inclined from the sliding vec- 
tor by angle a. The film flux, according to the 
lubrication theory, Is 

<t> = -( 1 h 3 grad p + 1/2 Uh [31] 

In the void dominated region liquid content Is 
given by the volumetric fraction or the fill fac- 
tor for each of the three morphological states; 
namely, © A , ©b» ©c* The velocity of the Internal 
boundary Is orthogonal to Itself. Therefore, the 
"sweeping flux vector' 1 Is 

Q * Vh {1 - <0 A + © B + e c >} [32] 

The transvoid liquid transport flux Is 

Y - Uh (6 a + 1/2 ©b> [33] 

Conservation of the liquid phase requires 

<<fr - ¥ - Q> • n =0 [34] 

Substituting Eqs. [31] to [33] Into Eq. [34], 
one obtains 

-<12n) -1 h 2 [(cos a) Op/ax) + (sin <x)Op/3z)] 

+ IKcos a) [ 1 /2 - <e A + 1/2 e B >] 

- v{i - <e A + e B + e c )} = o [35] 

Rearranging, one finds 

vn - (e A + e B + e c >] = u cos a [| - <e A + ^e B )] 

- (12jj) - ^ h 2 ^cos a |^ + sin a [36] 



• ►X 


Fig. 2. - Vector diagram. 


This is the equivalent Rankl ne-Hugoniot condi- 
tion. It shows that the sweeping speed of the 
Internal boundary depends on the neighboring film 
flux and the adjacent morphological states. A 
more restrictive form of this equation was 
originally presented by Olsson (8). 


IMPLICATIONS OF THE EQUIVALENT RANKINE-HUGQNIOT 
CONDITION 

Equation [36] applies to both the breakup 
boundary, where cos a > 0, and to the reformation 
boundary, where cos a c 0. It must be satisfied 
by all valid solutions of incomplete films. 

One of the solutions of the 01 sson-Ranki ne- 
Hugonlot condition Is the breakup line of a sta- 
tionary, Incomplete film. This corresponds to the 
Stationary breakup point: 

grad p = 0 [37a] 

© A + 1/2 © B - 1/2 [37b] 

The first of these is well known, but the second 
one Is rarely cited. 

Another solution of the stationary incomplete 
film refers to a point on the reformation line or 
Stationary reformation point; 

- (12 t i)" 1 h 2 {(cos a)(9p/3x) + (sin a) (3p/3z)} 

+ U(cos a){l/2 - (0 A + 1/2 e B >} = 0 [38] 

This relation was used by Jakobsson and Floberg 
(1) to find the shape of an enclosed cavity of a 
submerged journal bearing. 

For a moving Internal boundary, the 01 sson- 
Ranki ne-Hugonlot condition serves as a link be- 
tween an elliptic field of the fully filled domain 
and a hyperbolic field of liquid phase transport. 

In this context one may observe that, because the 
normal flux Is a unique consequence of the Dirlch- 
let problem, one can not rule out the possibili- 
ties that it may be Inwardly directed, and that 
the film pressure may be less than the boundary 
value. 

UNRESOLVED ISSUES 

In order to use the 01 sson-Rankine-Hugoniot 
condition, the three morphological state parame- 
ters © A , Ob, and ©c must be given. The ©c 
probably can be neglected unless there are special 
chemical causes for a biased affinity of the lubri- 
cant to the stationary bearing surface. Another 
plausible hypothesis Is to impose Eq. [37b] at 
the breakup point, irrespective of whether the 
internal boundary may be in motion. One is thus 
left with the remaining question of how liquid 
transport beyond the meniscus may be divided 
between an adhered film © A and streamers ©g . 
Suppose that they are given Initial values at 
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point E; then subsequent liquid-phase transport 
due to sliding can be stated as 

O/3t)[(0A - 0 B )h] + UO/a x)[(0 A + 1/2 0 B )h] - 0 

[391 

And, If one Is willing to neglect "cross flowing'* 
between the two constituents, the streamer trans- 
port is given by 

0/at)<e B h) + 1/2 UO/9x)<0 B h) = 0 [39a] 

and the adhered film transport Is given by 

(9/9tX0 A h) + U (3/3x) <e A h) - 0 [39b] 


Clearly, a physically permissible initial morpho- 
logical condition at E Is In the range 

0 S <0 A >e ^ 1/2 C40] 

If one wishes to Improve the speculative aspects 
of these Ideas, details of the three-dimensional, 
viscous, free surface flow with a moving surface 
and a moving attachment line must be studied. 

The Implementation of the Olsson-Rankine- 
Hugonlot condition In a numerical computation 
scheme is not trivial. It refers to the film flux 
at the boundary, which in turn shifts with time. 

It Is quite clear that the computation scheme must 
be capable of dealing with oblique mesh lines. 
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GEOMETRY AND MOTION OF THE MENISCUS 


JEAN FRENE, 4 LUC LAUNAY, 4 AND CODA H.T. PAN 5 

Dynamic loading is the primary concern in the 
design for journal bearings of automotive engines. 
Oil supply is to be optimized to ensure that the 
minimum gap Is adequate while friction is to be 
kept as low as possible. Operation with an incom- 
plete film is typical in such bearings. 

In the past the ir-film model is commonly 
used in analytical studies. As engineers strive 
for more accurate modeling tools, a more precise 
treatment of the dynamics of the incomplete film 
is In order. The ir-film model, as applied to the 
steady-state, plain journal bearing, is equivalent 
to the Gumbel approximation in rejecting the subam- 
bient portion of the pressure field. In conjunc- 
tion with the short bearing approximation, the 
ir-film model yields simple, closed-form analytical 
results for both static (eccentricity) and dynamic 
(radial velocity) effects. This combination of 
analytical tools has been the backbone for the 
study of journal bearings in automotive engines. 

The Reynol ds-Swi f t-Sti eber cavitation model (RSS) 
has been regarded to be more precise. Because it 
requires a null condition of the pressure gradient 
at the cavitation boundary, it had been thought to 
be intrinsically contradictory to the short bear- 
ing approximation and had been used primarily in 
"FDM type" numerical computations. More recently, 
a modified Reynol ds-Swi f t-Sti eber model was used 
to improve the short bearing approximation. Cavi- 
ties of both vented and enclosed types have been 
successful 1 y modeled . 

The Reynol ds-Swi f t-Sti eber model, which is 
restricted to the steady-state condition only. Is 
thus, valid for either a static load or a rotating 
load. As it was applied with the short bearing 
approximation, the results are appropriate for 
statid-load considerations of the journal bearing 
or for high rotating loads of the squeeze-film 
damper. The bearings of the automotive engine 
experience oscillatory loads which are not of the 
pure rotating type, and the static load component 
is not totally negligible. These applications 
motivate us to consider the more complicated prob- 
lem of a truly time-dependent Incomplete film. 

The short bearing approximation was retained for 
this work, so that one may direct full attention 
to issues related to the moving cavitation 
boundary. 

ANALYTICAL APPROACH 


The time-dependent, incomplete film problem 
requires the simultaneous treatment of an ellipti- 
cal domain and a hyperbolic domain, which are 
joined at the moving internal boundary where the 
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01 sson-Ranklne-Hugonlot condition Is to be satis- 
fied. A centrally fed short journal bearing will 
serve as the specific model in the present work. 
The dimensionless film thickness of a dynamically 
loaded journal bearing, neglecting angular mis- 
alignment, is 


H * 1 + e cos(9 - [41] 


In Eq. [41] both e and cj>, may be time-dependent, 
employing the short bearing approximation of 
Ocvlrk, the Poiseullle component of the circumfer- 
ential flux is neglected; thus, 

Y 0 * 1/2 coRCH [42] 

<|/ z - -(12^)- 1 (CH)3(a p /3z) [43] 


The continuity condition, with allowance for mo- 
tion of the journal , 1 s 

(3t 0 /9e> + R(8y /3z) + RCOH/at) = 0 [44] 


Upon making use of Eqs. [41] and [42], one finds 


(3 2 p/3z 2 ) » U2pC- 2 )H- 3 


Iff COS e - 4> 


b ■ 3t~ sin 9 ■ 


[45] 


Boundary conditions for the chosen problem are of 
the Dirichlet type. At the lubricant supply plane 

p(z = 0,9) * p s [46] 

and at the end of the film 

p<Z e ,6) * 0 [47] 

The end of the film may be at either the discharge 
end of the bearing or at the Internal boundary of 
the incomplete film; 1 ,e . , 

z e * L/2 [48] 

The elliptic problem can thus be determined by 
Integrating Eq. [45] with the constants of Integra- 
tion chosen to satisfy Eqs. [46] and [47], to yield 

p(z,0) = p <1 -2z/L> 


- 1/2 z(L/2 - z)<12p/C 2 )H 3 [e cos(0 - 4>) 

- (1/2 u - <j>) slnO - 4>>] [49] 
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LIQUID TRANSPORT IN THE DISCHARGE CAVITY 

In the discharge cavity liquid phase is 
assumed to be in the state of an adhered film. 
The liquid content is given by the adhered-film 
fill factor ffr. Transportation of the adhered 
film Is described by the equation 


But in the reformation branch f/\ is computed 
according to Eq. [51], so that Eq. [52] would be 
used In its full form. It is convenient to repre- 
sent the sweeping motion by 

(!t) z e " csc « V 


(3/at + to a/a©> Ha - o csoi 

where Ha =* fA H - The time-stepped representation 
of this solution 1 s 


- (1 - f A )"’[l/2 U <1 - 2f A > ctn a 

- (CK) 2 (12 v r'(£lp-|§ . S l» . |f)] £58] 


H A <0 f z;t) = H a (9 - «At, z; t - At) 

= 1/2 H(9 - <ot + cot, z; x) [51] 

for t > t, where t Is the instant the adhered 
film was formed at the breakup point. 

MOTION OF THE INTERNAL BOUNDARY 

The condition of 01 sson-Ranklne-Hugoniot 
gives the sweeping speed at an internal boundary 
to be 

V = (1 - f A ) -1 1/2 U(1 - 2f A ) COS a 

- (CH) 2 <12 H )" 1 ( ££ p-§ + sin a ||) C52] 

The sweeping speed V Is aimed along the cavity, 
pointing normal to the periphery. The boundary at 
any Instant is described by 

z e (©,t) ^ L/2 [53] 

which has a waist at 0 W such that 


Then updating of the boundary geometry can be 
executed by computing 

t 

z e <0,t> = z e (0, t - At) + (3/at>z e dt [59] 

J t-At 

at all inner points of the internal boundary. The 
end points move along z * L/2 and can be calcula- 
ted from 

<a/at)x = sec a V 

e z=L/2 z=L/2 

= (1 - f A > _1 1/2 IK 1 - 2f A ) 

- <CH> 2 <l2„>-'( R 2Eg-. tan «f§) ^ 

[60] 

Again, fA would be set to 1/2 for the breakup 
point but would be computed according to Eq. [50] 
from the prior history of the reformation point. 


<a/ae)zl Q = o [54] 

e'9, 


The inclination of the periphery normal is 

, -1/9 „ \ /0 [55] 

a * tan VR 36 z e) ~ * 2 

Branches on either side of the waist are 
joined with null slope. But analytic continuation 
from one to the other across the waist point can 
not be assured because computation of Eq. [51] 
would follow different rules for the two branches. 
In the breakup (0 < 0 W ), the Reynolds-Swift- 
Stieber law applies so that 

f A = 1/2 [56] 

and 

v breakup - -<CH) 2 <6p>- | (22p-| . sin . ff) 157] 


EQUATION OF MOTION 

The journal mass may experience an accelera- 
tion as the bearing film force may not cancel the 
applied force completely. This acceleration can 
be integrated In time to determine the varying 
velocity of the journal center. The correspond! ng 
squeeze-film terms would necessitate recalculation 
of the fluid film forces with the new cavity 
boundary. 

PERIODICITY 

In the automotive application the periodic 
state of the constant speed operation of the 
engine is of Interest. The character of the above 
procedure of calculation does not assure this peri- 
odicity in advance. Therefore, it is necessary to 
perform the computation for more than a full 
cycle, until the results demonstrate complete 
attenuation of the Initial transient caused by the 
assumed initial values of e and cj>. 
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CLASSIC CAVITATION MODELS FOR FINITE ELEMENT ANALYSIS 
J.F. BOOKER 8 


A review of several mathematical models com- 
monly used to incorporate cavitation In the finite- 
element analysis of lubrication is presented. The 
various mathematical models are completely equiva- 
lent In the limit, but they differ in physical 
Interpretation and numerical Implementation. The 
models are classic In the sense that they are not 
flow preserving. For clarity we review the devel- 
opment of finite-element lubrication analysis both 
with and without the complication of cavitation. 

CAVITATION EXCLUDED 

Transient lubrication analysis is usually 
based on a boundary-value problem in which the 
pressure satisfies the Reynolds differential equa- 
tion (30) over a complete film region, and either 
pressure value or normal pressure gradient, respec- 
tively, satisfy "essential and natural" conditions 
everywhere on its boundary, whose location Is fixed 
in the absence of cavitation. In an equivalent 
variational (extremum) problem (31), a pressure 
distribution Is sought which minimizes a func- 
tional (Integral) involving variants of the Rey- 
nolds operators over the region and "essential" 
boundary conditions on Its boundary. 

Following finite-element discretization of 
the pressure field, the functional can be expressed 
in terms of nodal pressure values. Necessary con- 
ditions for the resulting unconstrained minimiza- 
tion (quadratic programming) problem give a linear 
algebraic equation set. The set can also be ob- 
tained through application of Galerkin's method to 
the boundary value problem. 

Solution of the linear equation set Is com- 
monly accomplished by either direct decomposition 
(e.g., Gauss, Choleski, LUD (lower and upper 
decomposition)) or indirect iteration (e.g., 
Gauss-Seldel , and SOR) methods. 

CAVITATION INCLUDED 

In the classical view of transient cavitation 
analysis, the complete film region has one or more 
free boundaries on which both pressure value and 
normal pressure gradient satisfy so-called Rey- 
nolds boundary conditions (usually null). In a 
second view, a nonnegative solution (over the en- 
tire film region) is sought for the variational 
problem. The resulting constrained minimization 
(quadratic programming) problem no longer admits a 
simple linear algebraic solution. Third, we can 
choose to apply a mixture law view of two-phase 
(cavitatlng) fluids (15). Applied here, a com- 
plete film region is characterized by positive 
pressure and constant (mixture) density value 
(unity volume fraction); an incomplete (incipi- 
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ently cavitatlng) region is characteri zed by null 
pressure and negative density time rate of change. 
The physical Interpretation is obviously that of 
vaporization. It can be shown that the three view- 
points are equivalent In the limit. Any one is 
sufficient In principle to fix the (instantaneous) 
division of complete and incomplete film regions 
and/or the positions of their (free) boundaries. 

In practice, however, all three viewpoints intro- 
duce nonlinear complexities to the problem and a 
consequent layer of iteration to its numerical 
solution. 

Focusing on boundary location, as in the 
first view. Is particularly awkward for finite- 
element analysis. Since locations of free bounda- 
ries are unknown a priori, they must be selected In 
Iterative trials (with attendant nodal discretiza- 
tion) followed by solution of a well-posed problem 
over the complete film region (with only one of 
the two Reynolds boundary conditions fixed and the 
other evaluated for iterative closure). 

Focusing on local behavior Is particularly 
convenient for finite-element analysis, since it 
allows a fixed nodal discretization. The con- 
strained quadratic programming problem of the sec- 
ond view can be solved in a number of ways, all 
Iterative. Physically appealing is Chrlstopher- 
son's method, In which the linear equation set Is 
solved by Gauss-Seldel (or SOR) iteration, after 
which any computed negative pressures are set to 
zero, and the process repeated until satisfactory 
convergence is achieved. (One can imagine an 
analogous physical process in which quantities of 
liquid are converted to vapor.) The scheme has 
been used with success (32,33); its main drawback 
Is Its reliance on the iterative convergence of 
real variables. An alternative to the constrained 
quadratic programming problem is the linear comple- 
mentarily problem (equivalent in the limit). Its 
numerical application (34-36) appears identical to 
that of the mixture law (37), in which one seeks 
to partition nodes as cavitatlng (or not), depend- 
ing on whether particular equality and inequality 
constraints are satisfied. 

A metaphor would be a "porcupine bearing" in 
which all Interior nodes communicate via feed/bleed 
check-valves to an ambient pressure source. Al- 
though no cavitation at all would occur in such a 
preposterous arrangement, the mathematical problem 
is Identical to those above. 

Solution of all such problems is a matter of 
Integer programming (guesswork) combined with di- 
rect solution of linear equation sets for pres- 
sures and direct evaluation for flow terms. In 
Murty's method (analogous to Chri stopherson 1 s 
method) nodes at which constraint violations occur 
are shifted to the opposite partition for the next 
Iteration. The great advantage of such Integer 
iteration is that convergence is abrupt and unam- 
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biguous; with a bit of luck, it Is also fast. Fig- 
ure 3 shows pressure distributions from a series 
of such steps, starting with the assumption of no 
cavltated nodes whatsoever. It should be noted 
that such a blind start is unnecessary In practice; 
reasonable first Iterates result in convergence in 
one or two steps. (Corresponding nodal flow quan- 
tities do not provide for a meaningful plot in 
this previously unpublished study by LaBouff.) 

SUMMARY 


Several mathematical models commonly used to 
incorporate cavitation in the finite-element analy- 
sis of lubrication and suggested both physical in- 
terpretation and numerical implementation. These 
classic models, which are equivalent in the limit 
of discretization, describe incipient (quasis- 
tatic) ,nonevol vi ng (flow-preserving) transient cav- 
itation, in the sense that vapor only begins to 
form before the liquid is replenished. Though 
they are physically unrealistic in a really most 
basic way, the fact that they fix the cavitation 
situation quasi stati cal ly in instantaneous 
response to bearing kinematics means that the de- 
grees of freedom of the problem are hundreds of 
times fewer than if nodal density (volume frac- 
tion) must be called on to reflect evolutionary 
history. The models thus lead to computationally 
efficient algorithms and are, therefore, likely to 
persist for a long time to come, whatever their de- 
partures from realism. 


MATHEMATICAL MODELS AND RELATED ALGORITHMS IN CAVITATION PROBLEMS 

G. BAYADA, 7 8 M. CHAMBAT, 9 M. EL ALA0UI 9 


We concern ourselves with gaseous cavitation 
operating under both the steady-state and transient 
conditions. Applications to cl rcumferential ly sup- 
plied bearings and to seals are given. We describe 
first the widely used cavitation algorithms and the 
related mathematical problems especially In view of 
the output flow computation. We propose new algo- 
rithms using the finite elements variational ap- 
proach whose reliability has been theoretically 
proved. For unsteady problems, a time discretiza- 
tion procedure using the characteristic lines 
allows us to Introduce a set of stationary like 
problems for which equality between input and out- 
put flows Is preserved. 

STEADY STATE PROBLEMS 

Classical Models 

The Gumbel (or Half-Sommerfeld) model of 
cavitation is widely used because It Is easy to 
implement, but It Is associated to a false bound- 
ary condition on the free boundary surface. The 
situation is much more confusing for the Swlft- 
Stieber or Reynolds models. The validity of this 
last model has been supported in various ways: 

a. The regularity condition is required to 
hold at the free rupture boundary (38): 

p = 3p/3n * 0 C61 ] 

b. In numerous papers condition [61] is gener- 
alized to any free boundary (39). 

c. In a more mathematical formulation, condi- 

tion (61) is deduced from a minimization argument 
(40): Through all positive smooth functions satis- 

fying the external boundary conditions, the pres- 
sure is the function that minimizes the energy 

i ntegral . 

d. The pressure Is the solution of a varia- 
tional inequality (41,42). 

e. The pressure Is the solution of a comple- 
mentary probl em (43) . 

f. The pressure is obtained by the Christoph- 
erson procedure. 

There are two types of unpleasant problems 
that appear both analytically or numerically: 
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1. For a journal bearing with an axial supply 
line located In the convergent part of the gap, a 
strange pressure distribution has been obtained 
numerically with nonconnected active areas (44). 

2. The balance between the input flow and out- 
put flow is not always satisfied, especially for 
seals (45). 

Theoretical studies lead to the following conclu- 
sions (46): 

Formulations c, d, e, f are equivalent. They 
are associated with a well-posed problem 
which has a unique solution that fulfills the 
conditions of Eq. [61]. This solution is 
such that the cavitation area is located 
only In the divergent part of the mechanism. 

A New Model 

Various Interpretations of the model are 
described In references (1) and (14). A summary 
of those Interpretations Is given here. The pres- 
sure Is not a unique parameter governing the cavi- 
tation phenomenon; and we must Introduce another 
parameter 9 which can be viewed as a saturation 
function or a fractional film content (the percent- 
age of oil contained In an elementary volume) or 
the result of the homogenization of the oil fin- 
gers. In any case, the cavitated area is not 
empty as it contains some lubricant. 

The flow is now defined by 

<t> > h 3 /(12p)Vp - 6U h/2 [62] 

assuming (p ^ 0 and 0 < 9 < 1), with the con- 
straint that (9 - 1 If p > 0) then the condition 
d 1 v ($) » 0, In a weak variational sense, gives 
the equations for the problem. 

Numerical Method and Applications 

Elrod proposed an iterative procedure (15) 
well adapted for finite difference discretization. 
We propose a new procedure using the finite-element 
discretization, whose accuracy and convergence have 
been rigorously proved (46). 

It is possible to express the problem in terms of 
one single unknown X which is a solution of a 
nonlinear equation 

X = B(X) [63] 

This equation can be solved using the fixed-point 
method starting from zero. Once X Is known. It 
Is easy to recover the Initial functions p and 
9 . 

We consider a bearing with a circumferential 
supply line. We compute the load, the attitude 
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angle, and the Input/output flow rates for various 
supply pressures, various eccentricities, and 
three length to diameter ratios. The first compu- 
tation Is related to the Reynolds model solved 
using the Chri stopherson numerical procedure, and 
the second one to the Elrod-Jacobson-Floberg model 
using the new procedure. Both computations in- 
volved finite elements. Results demonstrate lit- 
tle discrepancy for the load between the two 
models, but gave more discrepancy for the eccen- 
tricity angle. The equality between output flow 
and Input flow Is very well preserved for the sec- 
ond model. However, it Is not true for the first 
model for which the discrepancy increases with 
eccentricity and length to diameter ratio. 

TRANSIENT PROBLEMS 

Classical Models 


The same analysis can be made for steady- 
state problems from both the theoretical and 
numerical aspects. The Chri stopherson procedure 
Is commonly used with an explicit or implicit time 
descretlzation. Bad results are usually reported 
for the balance of the flow, especially for 
seals. These will be Illustrated as a prime 
example . 


3(8h) 3_/V_ 3fi\ 1 3(8hU) 

3t 8h\ 12y 3x ) + 2 3x 


[ 64 ] 


This equation Includes all the above equations, 
especially the condition on the free boundary, 
which appears exactly as the classical Rankine 
Hugonlot condition for shocks in hyperbolic equa- 
tions and Is equivalent to those proposed by Elrod 
(14). Supplementary data involve (1) boundary 
conditions on the pressure and (2) initial data 
for 0. 

Numerical Procedure (47) 


The starting point is to substitute the total 
derivative 0(6h)/Dt for 1 /2 [ 9(©hu ) / 9x ] + (38h/at) 
with the meaning that the derivative in t is taken 
along the direction U so that the equation 
becomes 


D(8h) 3 /'h 3 3p^ n 

Dt 3x^l2y 3 xj ~ 


0=1 if p > 0 1651 


and then discretize the first term of the equation 
in time along the characteristic line. 

CONCLUSIONS 


New Model 

He describe It first In the particular case 
of an Infinitely long mechanism. However, the gen- 
eralization to real two-dimensional geometry is 
straightforward. 


An application to a seal was made in (47). 
Rectangular finite elements were used. 

1. The discrepancy between input mass flow 
and output mass flow never exceeds 4 percent In 
each period. 


As In the steady-state problem, a unique 
equation In a weak or distributional sense is well 
suited for numerical procedure: 


2. Periodicity of the load and of the flows 
Is obtained In four periods. The cavity follows 
the rotor movement when X 2 = 2xi . 
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NON-NEWTONIAN EFFECTS IN FILM CAVITATION 


A.Z. SZERI 1 0 , K. R . RAJAGOPAL 10 , AND ZHANG 10 

There exist In engineering usage a vast 
number of fluids* both naturally occurring and man 
made, the behavior of which cannot adequately be 
described by classical, linearly viscous theory. 
Though such fluids also abound In tribological 
practice, we have accumulated only scant knowledge 
of the effects non-Newtonian behavior might have 
on tribological performance. The principal prob- 
lem we are faced with In the area of non-Newtonian 
fluid mechanics Is our Inability to successfully 
model the stress-deformation relationship of a par- 
ticular material. But, In part, we can also blame 
the lack of reliable experimental data and, thus, 
the absence of guidance to theoretical research. 

Viscoelastic fluids differ from Newtonian 
fluids In one main aspect, viz., viscoelastic 
fluids possess memory. Having memory, for the 
fluid, means that the stress presently experienced 
by a fluid particle Is no longer determined by ex- 
isting conditions alone but depends on deformation 
history of the material In the neighborhood of 
that particle. The principal practical manifesta- 
tions of non-Newtonian behavior are normal stress 
effects and shear thinning. We have recently 
noted another strong effect on the velocity pro- 
file in a thrust bearing-like geometry (48) 

(Fig. 4). Velocity profiles were measured using 
laser Doppler anemometry, in a flow between rota- 
ting parallel disks. Figure 4 displays the tangen- 
tial velocity field at a fixed radial position, 
with the top disk stationary and the bottom disk 
rotated at constant angular velocity. The profile 
in (a) was obtained In water, and that in part <b) 
was measured in a viscoelastic fluid (0.029 per- 
cent polyacrylamide solution). It is easy to dis- 
cern the qualitative difference between the two 
profiles: Shear thinning occurs at a certain axi- 

al position, the separation distance of which from 
the bottom plate depends on various parameters, re- 
sulting in large velocity gradients. The cause 
for this localized action is not known to us. 

We have also noted recently (49-51) that jour- 
nal bearings lubricated with non-Newtonian lubri- 
cants could be less susceptible to thermal effects 
when compared with bearings lubricated with Newto- 
nian fluids. 

In this note we wish to further explore lubri- 
cation with viscoelastic fluids and focus our at- 
tention on cavitation, first wishing to place the 
constitutive relationship of the class of fluids 
we employ. In proper perspective. To this end we 
define the motion of a fluid particle as follows. 

If fluid particle X at time t occupies the spa- 
tial position x, we write 

x = X(X, t) t > 0 [66] 
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As t varies, so will vary the position x of 
the particle X. We insist that this relationship 
(Eq. [66]) be Invertible; l.e., that two particles 
do not at the same time occupy the same spatial 
position, and write for the Inverse of Eq. [66] 

X =■ x -1 <x.t) [67] 

This Is just mathematical shorthand for the state- 
ment that spatial position x at time t has the 




Fig. 4. - Azimuthal velocity between parallel, 
rotating disks, (a) Newtonian, (b) Viscoelastic 
fluid. 
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fluid particle X as Its temporary resident. 

Consider now another time instant x < t, at 
which time our fluid particle X was resident at 
the position £. Hence, we have 

£ = x<X,t> £68] 

But we can substitute for X from Eq. (67) into 
Eq. [681 and have 

£ * x[x _1 <^t) t x] E xt<x>t> [69] 

The relationship £ = xt (x * t) states that the fluid 
particle that Is In position x at the present 
time t, was In position £ at some former time 
Instant we call x. If we now let x ^ 0 vary 
while calling t the present time and keeping It 
fixed, we are looking backwards In time and trac- 
ing the history of motion of the fluid particle X. 

The vector valued function xt occurring In 
Eq. [683 Is called the relative deformation func- 
tion. If the deformation xt described is smooth 
and the material has no microstructure (and this 
applies to most deformations and most materials in 
engineering practice), we can define the gradient 
of xt » called F(s), through (52). 

d£ * F(s) dx; F(s) = 3£/3x; s - t - x [70] 

Equation [70] shows that the gradient of xt> 
our notation F(s), is Indicative of how an infini- 
tesimal line element d£ at the instant of time 
x had to be stretched and rotated so as to become 
the line element dx of the present time t. As 
we let s > 0 increase, the equality here signi- 
fying present time, we have the history of deforma- 
tion characterized by F(s). 

The deformation gradient history F(s) is 
used In the definition of a class of fluids called 
"simple fluids". Simple fluids (and the adjective 
"simple" refers to the lack of ml crostructure and 
does not imply that the fluid is elementary in any 
sense) are those In which the stress T at a 
point depends on the history of deformation of a 
neighborhood of that point, or 

s=«> 

T = f (F( s) ) [71 ] 

s=0 

There is a fairly good agreement among practition- 
ers of non-Newtonian fluid mechanics concerning 
the validity of Eq. [71] and the assumptions lead- 
ing up to it. But, there Is disagreement on how 
to proceed further, for Eq. [70] is hardly a 
practically useful statement. 

There are many ways to approximate the func- 
tional f In Eq . [7 1 ] . Integral representations 
are well received by the rheology community In gen- 
eral but have to be as yet handled numerically In 
a practical application. Differential representa- 
tions, on the other hand, yield to numerical treat- 
ment but allow only for infinitesimal memory. 


It has been shown that under certain condi- 
tions Eq. [71] has an expansion in terms of cer- 
tain differential quantities. When this expansion 
Is terminated at a low order and when account Is 
made of the particular geometry common to most 
lubricant films, we essentially obtain that the 
stress has the following structure (cf. (50)) 

T - - pi + [y + |l(tr Aj )] A, + oA 2 [72] 

Here, p is the molecular viscosity, a and 3 
are non-Newtonian parameters, and A} and A 2 are 
the first and second R1 vl 1 n-Eri cksen tensors, 
respectively, defined by 

+ •* T 

Aj = (grad v) + (grad v) 

[73] 

~ ~ t~ 

A 2 * dAj/dt + Aj(grad v) + (grad v) A 1 

The material parameter a promotes normal stress 
effects while a 3 < 0 portrays shear thinning 
and a 3 > 0 will make the fluid shear thicken. 
This is easily seen when Eq. [72] Is written for 
simple shear flow, for then 

2 

x xy =» [p + 3(du/dy) ] ( du/dy) [74] 

Models of this type (Eq. [72]) have been used by 
several researchers (53-55). We have applied It 
to journal bearings (50,51) and demonstrated that 
even a slight departure from Newtonion behavior, 
as measured by the numerical values of a and 3> 
renders bearing performance relatively Insensitive 
to changes in lubricant temperature. 

In a continuation of previous effort, we 
looked at non-Newtonian effects on cavitation in a 
journal bearing. Figure 5 shows the slow, gradual 
change in the pressure profile as the material 
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Fig. 5. - Pressure profiles in shear thickening 
fluid (y * -0.1). 
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Fig. 6. - Pressure profiles in shear thinning 
fluid (y = -0.1). 



Fig. 7. - Foldlines (S) on solution manifolds (M). 


parameter 3 Is made progressively smaller while 
a Is kept constant. To be precise, It Is not a 
and 3 that are used as the parameters of the prob- 
lem but their nondlmenslonal counterparts y and 
8, defined by 

y ; — 6 = [75] 

\xR \xC 

Here, U Is the surface velocity of a shaft of 
radius R, and C Is the radial clearance In the 
bearing. 

Figure 5 shows pressure profiles for shear 
thickening fluid, 6 > 0, at y = -0.1. As the 
value of S decreases from S * 0.004, the maxi- 
mum pressure, gradually decreases, otherwise the 
profile seems to change little. 

Figure 6 shows pressure profiles again at 
y a -0.1 but, now, for a shear thinning fluid. 
Starting at 6 =* -0.004, the pressure curve shows 
no cavitation. Increasing the value of 8 slowly 
Increases the maximum pressure, but, otherwise, 
the shape remains Invariant, still showing no cavi- 
tation. This phenomenon persists up to the value 
8 - 0.0025, a slightly shear thickening fluid. 

But at a somewhat larger value of S * 0.005, 

Fig. 5 already shows cavitation. 
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Fig. 8. - Pressure profiles at y * -0.1. 



Fig. 9. - Pressure maxima at y * -0.1. 


Our investigations are incomplete at this 
stage, and conclusions are not well founded. But 
they suggest the presence of a fold line on the 
solution manifold at 0.0025 < 8 < 0.005 when 
y » -0.1. Here, the two solutions, one with cavi- 
tation and one without, join. Of course it is pos- 
sible, or even likely, that the noncavi tat 1 ng solu- 
tion, If It at all exists, represents an unstable 
branch and Is therefore not reproducible in the 
laboratory. All we suggest at this juncture is, 
that there might be an other than the classical, 
cavitatlng solution to the problem. Our next task 
Is to trace the fold line on the solution manifold 
using continuation methods. Figure 7 Illustrates 
some possible fold lines of an equilibrium surface. 

Figure 8 shows two pressure profiles, one 
obtained with 3 « 0.0025 exhibiting noncavita- 
tion; while the other, calculated with 3 - 0.005, 
exhibiting a classical, cavitatlng profile. 

Except for the presence or absence of cavitation, 
the profiles are virtually Identical. Figure 9 
displays the pressure maxima at y » -0.1 as the 
parameter 3 is varied. 
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It Is well known that nonlinear problems 
often yield a rich multiplicity of solutions, sup- 
ported on Identical steady-state boundary condi- 
tions. As an example we may quote the work of 
Bejamln and Mullln (56), who found 21 solutions In 
the flow between short, concentric rotating cylin- 
ders. Their work was experimental and thus shows 


much more than Idle mathematical curiosity. This 
of course, constitutes no proof that what we have 
here In the non-Newtonian journal bearing Is also 
multiplicity of solutions; but It perhaps suggests 
that our findings should not be discarded out of 
hand. 
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CAVITATION MODELING FOR ANALYSIS OF FLUID FILM JOURNAL BEARINGS 

L.E. BARRETT 1 11 AND L. A. BRANAGAN 11 


A reasonable description of the cavitation 
phenomenon In thin fluid films is critical to the 
accurate prediction of the pressure distribution, 
and hence forces, in fluid-film journal bearings. 
Cavitation occurs when the fluid film is subjected 
to conditions of tension which tend to produce neg- 
ative pressures within the film. Usually, under 
these circumstances, the pressure falls to essen- 
tially the vapor pressure of the lubricant, but no 
lower. Gaseous components of the lubricant come 
out of solution and, together with vaporized lubri- 
cant, form a two-phase flow field of bubbles and 
streamers. Cavitation usually occurs in regions 
of diverging film thickness within the bearing and 
depends on the supply pressure of the lubricant to 
the bearing. Hence, cavitation regions can be 
drastically altered by changes In bearing geometry 
and supply pressure, and for a given bearing they 
are dependent on operating parameters such as 
shaft speed and bearing load. Figure 10 shows a 
comparison by Tonnesen of predicted and measured 
pressures obtained In a two-axial-groove journal 
bearing. This figure clearly shows the cavltated 
region of low pressure that occurs In the top pad 
of the bearing In which a diverging film thickness 
occurs. It should be noted that cavitation within 
bearing fluid films Is not considered to be Inher- 
ently bad. Only rarely has cavitation been cited 
as the primary cause of bearing damage or failure 
due to surface pitting, and cavitation control Is 
rarely a goal of bearing design except In the case 
of squeeze-film bearings. 

The fluid mechanics and chemistry of cavita- 
tion Is very complex. Fortunately, for most 
practical design considerations, a detailed 
description of the phenomenon Is not required. 

What Is required Is an accurate prediction of the 
cavitation region boundaries and the pressure, tem- 
perature, and shear heat generation within the 
region. This Is usually accomplished by applying 
reasonable auxiliary conditions on the fluid film 
pressures and temperatures obtained from a coupled 
Iterative solution of the governing equations, the 
Reynolds equation for pressure, and the energy 
equation for temperature. The solutions are depen- 
dent on boundary conditions of supply pressure and 
temperatures and may Include the heat conduction 
effects In the bearing and shaft. The simplest so- 
lution Is to solve the governing equations indepen- 
dent of the existence of cavitation and to set any 
predicted negative pressures to the cavitation 
pressure using a forward solution for the pressure 
starting at the leading edge of the film. It Is 
easily shown that this approach violates mass flow 
continuity, but it does provide a reasonable ap- 
proach to defining the cavitation region for many 
bearing applications, particularly if certain slm- 
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plifylng approximations to Reynolds equation are 
used, such as the short bearing approximation. It 
Is called the half-Sommerfel d pressure boundary 
condition . 

For finite-length solutions, this approach 
has usually been used to avoid an iterative solu- 
tion for the trailing boundary of cavltated region 
where the pressure Is the cavitation pressure and 
the pressure gradient is zero, thus ensuring mass 
flow conservation. This is referred to as the 
Reynolds boundary condition. Actually, the trail- 
ing cavitation boundary can be easily found with- 
out Iteration by back solving for the pressure, 
starting from the trailing edge of the fluid film 
and substituting the cavitation pressure whenever 
the solution predicts a subcavitation pressure. 

If coupled pressure-temperature solutions are 
being sought, the problem is already highly Itera- 
tive, so there is no advantage to be gained by 
using the half-Sommerfel d cavitation boundary con- 
dition, even if a forward substitution scheme Is 
employed. 

Once the cavitation boundaries are known, the 
heat generation within the cavltated region must 
be calculated. The heat generated has an indirect 
effect on the cavitation boundaries and no effect 
on the pressure within the cavltated region. The 
effect on the boundaries comes about because the 
heat generation alters the temperature within the 
fluid film in the cavltated region. Thus, the tem- 
perature of the lubricant exiting the trailing 
edge of the film Is altered. In turn, it affects 
the temperature of the film at the leading edge of 
the following pad. 

The usual theoretical approach is to use what 
Is called hot oil carry over calculations. In 
this method the lubricant exiting the trailing 



Fig. 10. - Predicted and measured pressures in a 
two-axi al -groove bearing. 
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Fig. 11. - Predicted and measured pressures in a 
two-axial -groove bearing. Experiemnt 
conditions: Tonnesen bearing; 3500 rpm; 5600 N; 

T s * 50 °C. Present work used 2° elements and 30 
harmonics . 



TOP VIEW 

• CAVITATION ONSET 



SIDE VIEW 



Fig. 13 . - Predicted and measured temperatures in 
a two axial groove bearing with two-phase flow 
model used in cavitated region. Effect of 
cross-film bearing; Tonnesen bearing; 3500 rpm; 
Babbitt temperature profile; 1289.5 Ibf load; w / 
thermal cav. 

loped in the cavitated region. This model treats 
the cavitated region as one with a central streamer 
of lubricant flowing along the bearing centerline. 
The heat generation is easily calculated since the 
velocity gradient across the film is essentially 
constant, at least for laminar flow, being approxi- 
mately Rw/2, where R is the shaft radius and 
cu is the shaft angular velocity. 

Calculations performed using this model are 
shown in Fig. 11 for a two-axial-groove bearing. 
This figure shows the circumferential temperature 
variation of the bearing surface in both the 
loaded pad and cavitated pad. Measured tempera- 
tures made by Tonnesen are shown. Although heat 
conduction into the shaft and bearing surface is 
included in the theoretical calculations, the 
analysis predicts significantly higher tempera- 
tures in the cavitated region than those meas- 
ured. This suggests that the heat generation 
model is predicting larger shear heating than 
actual ly occurs. 


Fig. 12. - Model for two-phase flow in cavitated 
1 ubricant films. 

edge of a pad is assumed to cling to the shaft and 
to be carried into the inlet of the following 
pad. If this volume of lubricant is Insufficient 
to fully supply the following pad, the remaining 
volume is made up of lubricant supplied through 
the Inlet supply groove. The temperature of the 
Inlet lubricant is determined by the ratios of the 
carry over and the supply groove volumes at their 
respective temperatures. 

A detailed theoretical model of the two-phase 
flow in the cavitated region is usually not con- 
ducted. Some assumption of the flow field is 
therefore made, and the shear heating of the lubri- 
cant In this region is determined by supplemental 
calculations. The simplest model is to assume 
that the lubricant film extends from the bearing 
surface to the pad surface. The axial film width 
decreases in the diverging circumferential direc- 
tion since there is a constant volume of lubricant 
and because there is no pressure gradient deve- 


Figure 12 shows a two-phase lubricant-gas 
model for heat transfer in the cavitated region In 
which the lubricant adherence to the shaft surface 
can be assumed. Further, it is assumed that the 
height of the film diminishes along the diverging 
portion of the gap rather than have the film 
decrease in width. This model has been developed 
by observations of flow fields in transparent bear- 
ings, which suggest that this type of flow pattern 
exists. Because the viscosity of the gaseous por- 
tion of the film is considerably lower than that 
of the lubricant, most of the shearing occurs in 
the gaseous region. Thus, this model generates 
less shear heating in the cavitated region. The 
heat conduction to the shaft and bearing surfaces 
is also modified because of the difference in 
thermal properties of the two fields. 

Figure 13 shows the results of calculations 
performed on the same two-axial-groove bearing as 
presented in Fig. 11. The results show better cor- 
relation with the measured temperatures in the cav- 
itated region, suggesting that the two-phase model 
used Is reasonably good. 


48 


A MODEL FOR CAVITATION IN DYNAMIC BEARINGS 

JOHN A. TICHY 12 AND C.P. KU 12 


In fluid film lubrication problems, pressure 
often seeks a level lower than lubricant vapor 
pressure In the divergent region of the bearing, 

In which case gas or vapor forms a cavity. For 
the steadily loaded bearing, the cavity Is usually 
formed due to ventilation from the surrounding 
atmosphere or the emission of dissolved gases from 
the lubricant. For the dynamically loaded bear- 
ing, the cavity Is composed mostly of lubricant 
vapor. 

For the heavily and moderately loaded bear- 
ing, the Jakobsson-F1oberg-01 sson cavitation model 
seems to be correct to predict the extent of the 
cavitation region and to specify proper pressure 
boundary conditions on the Interface. However, 
this model Is difficult to apply due to the nonli- 
nearity of the reformation boundary conditions. 

It requires extreme numerical effort In guessing 
the shape of the cavitation region and checking 
the flow balance (57) . 

In the proposed model we assume the shape of 
the pressure axial profile, and the Floberg bound- 
ary conditions are maintained. After integrating 
the Reynolds equation across the axial direction, 
a nonlinear ordinary differential equation for the 
distance of the cavitation boundary from the bear- 
ing centerline as a function of the circumferen- 
tial coordinate Is obtained. In the present study 
we consider the circular-centered orbit squeeze- 
film damper bearing with two circumferential oil 
supply grooves at the bearing ends (Fig. 14). The 
lubricant properties and temperature are assumed 
to be constant, and the fluid Is laminar and Iner- 
tialess. The unwrapped dimensionless bearing sur- 
face and pressure boundary conditions are shown In 
Fig. 15, where Z * 0 Is the centerline and W(X) 
Is the width of the cavity. 

ANALYSIS 


The dimensionless Reynolds equation and the 
Floberg boundary conditions are 
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where C Is shown In Fig. 15 X - L'/Rq, and X] , 
x 2» x m> W(X) are unknowns. The P(X,Z) Is a peri- 
odic function In the X direction and global mass 
conservation requires 
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[78] 


After Integrating the Reynolds equation 
across the Z direction and applying boundary con- 
ditions at Z - 0 and Z - 1 or W(X), an Integral 
Reynolds equation In regions I, II, and III Is 
obtained: 

Region I: 



10 dH 12,, 
12 dX + ~>?Zl 


[793 



Fig. 14. - Squeeze-film damper schematic. 



Fig. 15. - Bearing surface. 
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Region II: 

|L(h 3 f) . 12(1 - »>f . % ZL [SOI 
Region III: 

a*(" 3 l)- ' 2<l - ">35*%.* ,2IH< «> - HI ® 

[81] 

where the Integral pressure P(X) Is defined as 
the Integral of P(X,Z) across the Z direction. 

Equations [82] and [83] are pressure profile 
assumptions In the noncavitation and rupture re- 
gions, respectively. These pressure profiles sat- 
isfy all boundary conditions at Z = 0,1 and Z - 
H. 


P N <X,Z) - P S Z 2 +■ P C <X><1 - Z 2 ) [82] 



[83] 
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Fig. 16. - Cavity shape. 
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|~<X,H) + ||(X,W)^ - o [84] 

Calculating the Integral pressure of Eqs. [82] and 
[831 and applying the chain rule (Eq. [84]), we 
substitute these equations Into Eqs. [79] to [81] 
and apply the boundary conditions at X], X m , and 
X 2 on the Interface to obtain the following ordi- 
nary differential equations: 

Region I: 


3 dH^c 
H dX dX 


X 2 (P C ' V 


W ■ dT <X l> - P c (X 2 ) 


Fig. 17. - Centerline pressure. 
Region III: 
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Region II: 
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6 1 36(1 - H) dH n 

X 2 1 - W + P $ H 3 dX ' 


W(X, > = ^(X ) > 0, ^<x.) 
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In Eqs. [85] to [913, the independent variable Is 
X; the dependent variables In the three regions 
are Pq, W, and P; and H Isa known function. 
Known parameters are P$ and X. Unknown parameters 
to be determined from the solution procedure are 
the locations X], X m , X 2 (see Fig. 15) and the 
slope dW/dX at X] (= DW1 ) . 
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SOLUTION PROCEDURE 


RESULT 


First we guess the values of X] and DW1 . 
Equations C85] and [87] are solved with the ini- 
tial conditions at X - Xj , and the values of X 2 
and X m are determined with extra boundary condi- 
tions in Eqs . [86] and [88]. In principle, only 
one boundary equation Is needed for each inter- 
face. With the Initial conditions obtained from 
Eq. [87] at X - X m , we solve Eqs. [89] and [90]. 
The guessing values of X] and DW1 are adjusted 
until Eq. [91] and Eq. [78] are satisfied. 


The typical cavity shapes and the centerline 
pressure distributions are shown In Figs. 16 and 
17. Compared with the numerical schemes, with the 
current model we only need to guess two values 
rather than the whole cavity shape, and the 
partial dlfferentail equation reduces to odinary 
differential equations. A large amount of numeri- 
cal work is saved and not too much accuracy Is 
lost. We are currently obtaining the numerical re- 
sults by finite-element methods. Compared with 
the short bearing model (58), the present model ob- 
tains better results when correlated with the 
finite-element-method model. 
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THE INFLUENCE OF CAVITATION ON JOURNAL- BEARING ROTORDYNAMIC CHARACTERISTICS 


M.l. ADAMS 13 


Even under the smoothest running conditions, 
the journal always exhibits some orbital vibration 
about Its equilibrium position with respect to the 
bearing. Provided this vibration is sufficiently 
small, one can justify characterization of the 
resulting fluid-film Interactive dynamic force 
between journal and bearing by a linear mathemati- 
cal model. The essential Ingredients of this model 
are the linearized stiffness and damping coeffici- 
ents for the journal bearing, l.e., a 2 by 2 stiff- 
ness matrix and a 2 by 2 damping matrix. These 
rotordynamlc coefficients for the Interactive dy- 
namic force between the journal and bearing are de- 
lineated in the following equation: 


C92] 


The radial dynamic force f exerted on the jour- 
nal is the vector sum of the applied static force 
W and the total fluid-film force F, as Illus- 
trated In Fig. 18. 

In the simplest of cases, the flow field be- 
tween the rotating and nonrotating members has no 
memory. This would be the situation for the hydro- 
dynamic pressure distribution produced In a purely 
viscous, fully developed, Incompressible fluid- 
film journal bearing In the absence of cavitation 
(film rupture) as characterized by the following 
classical lubrication equation of Reynolds com- 
bined with pressure boundary conditions: 




where n Is the unwrapped circumferential coordi- 
nate and £ is the axial coordinate map of the cy- 
lindrical bearing surface. 

The lubricant fluid film pressure distribu- 
tion p(n»B,t) embodies an interactive radial 
force vector which acts to separate the journal 
and the bearing. This pressure distribution Is 
functionally related to the radial-orbit position 
and velocity vectors of the journal center rela- 
tive to the bearing center; that is, the film 
thickness distribution h separating the journal 
and the bearing Is a function of the relative posi- 
tion of the journal with respect to the bearing as 
follows : 


h = C - x cos(n/R> - y s1n( n /R> C94 J 


1 department of Mechanical and Aerospace 
Engineering, Case Western Reserve University, 
Cleveland, OH. 


where x - x s + x and y « y s + y (see Fig. 18). 
Equation [94] Is applicable to the typically 
assumed condition that bearing clearance is much 
smaller than bearing radius (C « R) and that the 
bearing and the journal axial centerlines are par- 
allel; i.e., no ml sal Ignment, as given by 3h/3£ 

* 0. Consistent with Eq. [94], the following 
derivatives of h, which appear in Eq. [93], are 






(b) 

Fig. 18. - (a) Radial motion coordinates of jour- 
nal relative to bearing, referenced to static 
equilibrium state (W = static load); (b) Harmon- 
ic orbital motion of journal relative to bearing; 
x := x 0 slnCQt + e x ); y = y 0 s i n (at + 9y). 
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written as 

dh/dn - [<x s1n( n /R) - y cos( n /R>]/R [95] 

dh/dt - -x cos(n/R) - y s1n(n/R) C 96 J 

Clearly, In this case, the Instantaneous pressure 
distribution and, thus, the Instantaneous journal- 
bearing separating force are determined by the 
Instantaneous position and velocity vectors (x, y) 
and (x, y) , respectively, and thus, has no memory 
feature. With no memory, the eight stiffness and 
damping coefficients are functions only of the 
static equilibrium flow field, e.g., for purely 
laminar flow, a function of bearing Sommerfeld 
number only. 

When cavitation Is admitted into the problem, 
the longstanding approach is to require p « p va p 
and 9p • n - 0 (to conserve mass flow) at the 
Interface boundary separating the full-film region, 
In which Eq. [931 Is solved, and the cavltated re- 
gion, In which all pressures are set equal to the 
vapor pressure, p V ap* this *PP r °ach to cavl- 

tatlng bearings, the solution again leads to rotor- 
dynamic coefficients that are functions of the 
static equilibrium flow field only. While this ap- 
proach to cavitating bearings Is probably an ac- 
ceptable characterization for static load capacity 
predictions, It has, unfortunately, been extended, 
with little questioning, to the rotor-dynamic coef- 


ficients. If closely scrutinized, the approach 
for cavitating bearings Is found to satisfy mass 
conservation only within the full-film region and 
on Its boundaries but not within the cavltated or 
ruptured region. In essence, this approach func- 
tions within a rotordynamlc context In the follow- 
ing loosely worded manner: whenever the right-hand 

side of the Reynolds equation gives the command, 
the left-hand side instantly becomes activated to 
allow the generation of local film pressures, 
whether or not the fluid has had time to refill 
the local gap. 

In fact, the actual system is quite transient 
or memory-like, as shown by the computations of 
Elrod and Adams (14), who employed a time-transient 
solution scheme which properly accounted for mass 
conservation throughout the bearing clearance gap. 
Thus, even In the absence of fluid Inertia, a dy- 
namically loaded cavitating journal bearing can ex- 
hibit Inertia-like effects. And these effects can 
be expected to become progressl vely more Important 
as vibration frequency is increased. This is 
quite similar to compressibility effects. 

A more complete discussion of this and related 
topics Is given In (59). For the purpose of this 
discussion, It Is simply concluded that the memory 
or frequency effects of cavitation on bearing 
rotordynamlc characteristics Is important. 
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CAVITATION IN JOURNAL BEARINGS AND SQUEEZE-FILM DAMPERS: ANALYTICAL 

B. JACOBSON 1 4 


Cavitation boundary conditions for stationary 
and vibrating journal bearings have been discussed 
In the literature for more than 30 years. The 
Swift-Stleber condition for the start of the cavi- 
tation region has been used both for stationary 
bearings and whirling bearings where the minimum 
film thickness is constant and the speed of the 
cavitation zone Is less than the mean value of the 
surface speeds. A more accurate boundary condi- 
tion for the start of the cavitation zone was 
given by Floberg where the number of streamers per 
unit width in the cavitation zone determined the 
amount of subcavity pressure In front of the cavi- 
tation zone. The most frequently used outlet 
boundary condition for the cavitation zone Is the 
JFO (Jakobsson, Floberg, Olsson) condition, which 
says that the oil flow going into the cavitation 
region Is transported straight through It without 
side flow and that there is continuity of mass 
flow at the end of the cavitation region. 


V = 


W.i (h - r) 

V h 


[98] 


where h c is the film coefficient of heat trans- 
fer for condensing oil vapor at low flow velocity, 
tv Is the heat of evaporation for oil, and AT v j 
is the temperature drop from vapor to liquid oil! 

For a 10-mm-dlameter , flat bubble In a 20-pm 
oil film and with p = 6 x 10*6 kg/m^, h q * 600 
W/m2.°C, / v = 350 kJ/kg, and the collapse velocity 
V * 5 m/s, the temperature drop needed from the va- 
por to the oil Is T v j * 7.0x10*5 C. This means 
that the vapor bubbles collapse almost without any 
resistance. 


CAVITATION FORMATION 


THE DIFFERENCE 8 ETWEEN VAPOR AND GAS CAVITATION 

Whereas the pressure In a vapor bubble is 
almost totally given by the temperature, the 
pressure In a gas bubble increases rapidly accord- 
ing to the gas law as the volume of the bubble 
decreases. The transport of gas into and away 
from a bubble Is a slow process, as can be clearly 
seen from the cavitation films where It often 
takes hundreds of vibrations before the gas cavita- 
tion pattern stabilizes. If the gas transport 
Into and out of the bubble Is neglected, the bub- 
ble pressure can easily be calculated using a poly- 
troplc compression, i.e., midway between an Iso- 
thermal and an adiabatic compression. This gives 
a radius ratio of 0.2 for a flat bubble when Its 
pressure Is increased 100 times. Therefore, If 
the bubble contains air, It will be visible in the 
high-pressure zone of the bearing. Vapor bubbles, 
on the other hand, collapse within milliseconds of 
formation. 

ANALYSIS OF VAPOR CAVITATION 


A simple model of the vapor cavitation behav- 
ior in dynamically loaded journal bearings is to 
assume the pressure to be constant during the col- 
lapse of the bubble. This means that the vapor 
condenses at the same rate as the bubble wall 
moves. For a spherical bubble, this gives the 
col lapse veloci ty as 

MT , 

V = [97] 

V 

and for a flat bubble with the radius r in an 
oil film of the thickness h, 


15 Chalmers University of Technology, Gote- 
borg, Sweden, and SKF Engineering and Research Cen- 
tre, Nleuwegein, The Netherlands. 


By changing the coordinate system for the 
vibrating bearing to a system rotating with the 
same speed as the minimum oil film thickness, the 
right-hand side of the Reynolds equation can be 
simplified for times close to the Instant the mini- 
mum oil film thickness has its lowest value. The 
right-hand side of the Reynolds equation for the 
geometry of the vibrating bearing is then simpli- 
fied to 


X = C 


sin 9 


V e s * E d ) 
2 


e d“d, 


[991 


For close to zero minimum oil film thickness, e s + 
ejj + 1 , and the criterion for the formation of 
vapor cavitation Is simply 


e d 



[ 100 ] 


Just as was mentioned In the experimental results, 
a high rotating speed u s of the journal stops 
the vapor cavitation (see also (24)). This Is 
easily seen as must always be less than 1 . 


CAVITATION COLLAPSE 


Once vapor cavitation has been formed and 
some time has passed, the minimum oil film thick- 
ness Increases, and the atmospheric pressure at 
the edges of the bearing becomes important for con- 
sideration of the collapse of the vapor cavitation 
bubbles. The side flow induced by the atmospheric 
pressure at the bearing edges increases with the 
third power of the oil film thickness. The high- 
speed films showed that the bubble collapse is 
mainly governed by this side flow. A simple short 
bearing type of analysis can give the criterion 
for when the bubble reaches Its maximum width and 
starts to collapse. This happens when 
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(3) ROTOR AT MAXIMUM SEPARATION (t = 0; 
£ = 0 . 10 ). 



1 - MAXIMUM FILM 

(d) BUBBLE EXPANDS (t = 42.51 MS; 
£ = 0.83). 




L MAXIMUM FILM 

(e) BEGIN SEPARATION (t = 33.89 MS; 
£ = 0.90). 



(C) CAVITATION INCEPTION U = 25.41 MS; 
£ = 0.83). 



^-MAXIMUM FILM 


(f) OCCURRENCE OF PEAK PRESSURE 
t = 29.31 MS; £ = 0.88) . 


Fig. 19. - Pressure distribution and bearing configuration for a full period 
of shaft whirl (0.1 i e i 0.9). Figures (a) to (f) viewed clockwise are con- 
secutive in time. 
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If the cavitation collapse takes place In the same 
location many times, it will result in cavitation 
erosion damage on the journal and the bearing. 


NUMERICAL-THEORETICAL CALCULATIONS OF 
VAPOR CAVITATION 


The appearance and motion of the cavitation 
zone and Its boundaries can be studied using numer- 
ical methods. This was done using a computer pro- 
gram (7) In which Elrod's algorithm for the conti- 
nuity of mass flow into, through, and out of the 
vapor cavitation zone was used. This program makes 


it possible to predict the dynamic behavior of the 
vapor cavitation zone. The Cray computer at NASA 
Lewis Research Center was used to do the computa- 
tions. The boundary condition for the vapor cavi- 
tation zone was the pressure equal to zero, which 
is consistent with the earlier experimental re- 
sults. This was used together with Elrod's algo- 
rithm to make sure that the continuity of mass 
flow was maintained, both in the full oil film and 
In the cavitation zone, during the whole vibration 
cycle. Some numerically calculated pressure dis- 
tributions and cavitation regions (23) are shown 
in Fig. 19. A comparison with the experimental re- 
sults shows similar behavior for the experimental 
cavitation zone and the theoretically calculated 
zone, even if the size of the experimental cavita- 
tion zone is slightly smaller. This can be partly 
explained by elastic distortion of the PMMA tube. 
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CONCLUSIONS 

Vapor cavitation occurs when the tensile 
stress applied to the oil exceeds the tensile 
strength of the oil or the binding of the oil to 
the surface. The physical situation necessary for 
vapor cavitation to exist Is a squeezing and slid- 
ing motion In a bearing. The Implication of the 
results as they apply to squeeze film dampers Is 
that If the sliding velocity Is Increased, the 
occurrence of vapor cavitation Is decreased (see 
also (24)). 


Vapor cavitation bubbles collapse at pressures 
lower than the atmospheric pressure, and they can- 
not be transported through a high pressure zone of 
the bearing, unlike gas cavitation bubbles. 

The numerical calculations (23) of the cavita- 
tion bubble dynamics show that an increasing vibra- 
tion ampl 1 tude gives: 

• Longer vapor bubble life 

• Larger maximum vapor bubble size 

• Earlier start and later collapse of the va- 
por bubble. 
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The following discussion summarizes the 
major issues brought out by the contributors as 
these issues relate to present technological 
requirements. At the risk of appearing myopic, 
the editors have taken the liberty to occasion- 
ally emphasize the importance of the work, what 
has been resolved, and what needs to be done. 

Many of the articles focused on cavitation 
as it pertained to a broad list of applications 
in crankshaft bearings, rotordynami cs , slider 
bearings, and seals. The key issue is the 
dynamic loading. As pointed out in the intro- 
duction to Part II, dynamic loading causes 
changes in the local film thickness, which leads 
to nonstationary cavitation. This certainly 
plays an important role in determining rotordy- 
namic characteristics. Adams' discussion ampli- 
fies this topic, especially concerning the 
effect of frequency. Adams draws the analogy 
that nonstationary cavitation can exhibit 
"inertia-like effects" and that these effects 
can be expected to become progressively more 
important as vibration frequency is increased. 

Of course, this is only one example of many 
where nonstationary cavitation is likely to have 
a significant effect on the rotordynamic charac- 
teristics. Jacobson discusses the effect of 
deformation, vibrational amplitude, surface ten- 
sion, and sliding velocity on vapor bubble 
inception, size, life, and collapse. Both 
Inception and collapse play an important role in 
the bubble dynamics which, in turn, affects 
rotordynamic behavior. 

As pointed out by Launay, dynamic loading 
is the primary concern in the design of journal 
bearings of automotive engines. Here, operation 
with an incomplete film is typical in such bear- 
ings, and one strives to maintain an adequate 
film while keeping the friction as low as pos- 
sible. In the case of diesel engines, the com- 
bustion conditions can become more severe (i.e., 
higher rates of change of cylinder pressure), 
resulting in cavitation erosion damage to the 
bearings. Martin discusses the work on this 
subject by Garner, James, and Warriner (52). He 
Includes examples of typical cavitation erosion 
damage, together with a plausible mechanism of 
damage. Commonly applied palliatives and their 
effectiveness, based on engine experience, are 
also discussed. It was concluded that cavita- 
tion erosion damage resulted from three main 
causes. Two were related to oil-column inertia 
effects in the connecting rod drillings, crank- 
shaft drillings, and bearing oil feed features. 
These were obviated by simple modifications to 
the detail design. The third, related to jour- 
nal movements within the bearing clearance 
space, posed a more difficult problem. Martin 
indicated that a more detailed theory than that 
used in their analysis was needed to reliably 
predict damage for the design stage. 

Of course, cavitation effects are not lim- 
ited to journal bearings and squeeze-film damp- 
ers. In a seal or slider application, Lebeck 




describes a mechanism in which cavitation is 
used to advantage. The asperities in parallel 
sliders were hypothesized to produced load sup- 
port by virtue of the fact that cavitation deve- 
loped at the downstream side of the asperity. 

In the presence of cavitation there is a net 
load support for individual asperities and, 
hence, for the parallel slider. 

The above articles dealt mainly with cavi- 
tation effects in various applications, especial 1 
where it concerned dynamic operating conditions. 
In order to gain insight and provide better 
design through numerical modeling for these 
applications, we must identify the modes of cav- 
itation associated with different dynamical con- 
ditions. For example, in the introduction to 
Part I we distinguished between gaseous and 
vaporous cavitation, but in actual practice 
both may occur together unless definitive steps 
are taken to preclude one or the other (i.e., 
degassing the fluid so that only vapor cavita- 
tion can occur or pressurizing the fluid so that 
vapor pressures are not possible). Jacobson 
elaborates on some of the noted differences 
between vapor and gaseous cavitation in Part II. 
Parkins' article on the "Categorization of Cavi- 
tation Types and Onset Conditions" is a begin- 
ning in understanding the conditions under which 
one or the other might dominate. The categori- 
zations in this article were based on sustained 
oscillatory squeeze motion - simplified by 
removing the translational component. Three 
cavitation regimes - including subregimes - were 
Identified and characteri zed according to the 
method of feeding the bubble. It was noted that 
cavitation bubbles fed from within the film 
could be either gaseous or vapor. Cavitation 
bubbles fed from outside the film obviously were 
gaseous (e.g., refer to the proposition of "side 
entrainment of gas," by Heshmat) and took sev- 
eral forms. Three of these forms were leaf type 
patterns, and the fourth appeared as smoothly 
curved bubbles entering from all sides of the 
square plate. Certain regimes were observed to 
coexist. Parkins indicates that "a limited 
knowledge is available about the conditions which 
determine the onset of a particular regime or 
subregime and its transformation into another 
subregime. The large number of possible forms 
of behavior make a confused picture when 
observed over a small range of variables, which 
renders it difficult to formulate classifica- 
tions." In the editors' opinions, this work 
deserves an expanded effort in which classifica- 
tion according to methods of feeding the bubble 
should be looked at with the translational com- 
ponent as well. It is only through the classi- 
fication procedure that we can hope to bring 
order to this "confused picture," thereby pro- 
viding a basis upon which analysts can formu- 
late appropriate models. 

We have already noted the importance of 
bubble dynamics in rotordynamic behavior. Bub- 
ble dynamics is a complicated subject in itself. 
The dynamics of the bubble can be affected by 
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the gas and vapor content (or fractional film 
content within the bubble) as well as by trans- 
cavity flow rate or fluid propagation speed 
within the bubble. The propagation speed of the 
fluid through the bubble depends on the morpho- 
logical condition of the bubble. Pan and Frene 
discuss this issue with regard to numerical 
modeling and implementing the transport process 
across the interface. They include the possi- 
bilities that a portion of the void-dominated 
region is locally filled with liquid, oridging 
both surfaces, or that a thin liquid layer may 
be attached to either but not both surfaces 
(also see discussions by Heshmat by Barrett and 
Branagan in Part I). For the first possibility, 
the propagation speed of the fluid is assumed 
to be the average of the two surface speeds. 

But for the latter possibility, the fluid takes 
on the speed of the surface to which it becomes 
attached. In the case of the journal bearing, 
one might argue that fluid transport through the 
cavity influences the position, size, and life 
of the bubble. Consequently, it can have an 
effect on the load capacity (7) and, possibly, 
the stiffness and damping character i sti cs of the 
bearing. Knowledge is scant on the conditions 
that determine whether attachment is to one or 
both surfaces, although it is generally accepted 
that load is a determining factor. More work 
is needed to clarify this issue and the effect 
it has on bubble dynamics and, consequently, on 
rotordynamic performance. 

The temperature effects on cavitation are 
an important consideration both from the stand- 
point of bubble dynamics and heat-transfer 
characteristics. Both Heshmat and Barrett and 
Branagan discussed temperature effects as it 
pertained to "hot-oil carry over" in padded 
bearings. Not surprisingly, it was noted by 
both that heat generation within the cavitation 
region was practically nonexistant. Barrett and 
Branagan compared theory with experiment. They 
found that treating the void region as a two- 
phase lubricant-gas in which the fluid attaches 
itself to the shaft was the most suitable model 
for the experiment considered. Further studies 
involving a wider range of operating conditions 
would be useful information to the analyst in 
modeling thermal effects on cavitation. Cor- 
rectly predicting the apportionment of heat 
transfer to the shaft and housing is partly 
dependent on choosing the appropriate cavita- 
tion model. One other way in which thermal 
effects can affect cavitation is through the 
indirect changes in surface energy. In this 
context we refer to temperature effects on the 
wettability of the lubricant to the surface and 
the consequences to fluid transport through the 
bubble as well as heat transfer to the bearing 
and shaft material. Albeit, Szeri , Rajagapol , 
and Zhang (Part II) noted that journal bearings 
lubricated with non-Newtonian lubricants could 
be less susceptible to thermal effects when com- 
pared with bearings lubricated with Newtonian 
fluids. 


Non-Newtonian effects on cavitation were 
further discussed by Szeri and shown to be 
important. That is, a shear thickening fluid 
cavitated within a range of material parameters, 
whereas the shear thinning fluid did not cavi- 
tate at all for the same parameter range. 
Furthermore, the change in material parameter 
for the shear thickening fluid affected the 
position of th rupture boundary to some extent. 
Since many of the lubricants used these aays are 
non-Newtonian, it behooves us to learn much 
more about this phenomena. 

Using high-speed photography, Jacobson and 
Hamrock (Part I) noted a clear influence of the 
surface energy of the journal surfaces on cavi- 
tation incipience. For steel journals, cavita- 
tion bubbles start in different parts of the 
low-pressure region, depending on where the 
largest bubble or dirt particle is located. For 
PTFE journals, the cavitation seems to start 
from the same position of the journal for each 
vibration, which indicates a rupture of the oil- 
PTFE interface. Presently, analytical models do 
not account for this phenomena. 

This calls attention to the subject of 
numerical modeling. Most of the contributions 
to this subject have dealt with flow preserving 
schemes, numerical expedience, and some recent 
advances using finite-element schemes. In flow 
preserving schemes, the attention is naturally 
drawn to the more complex cavitation boundaries. 
With the advent of the Elrod algorithm, much of 
the complexity has been reduced while retaining 
the flow-preserving features of the JFO theory. 
Keith and Vi jayaraghavan (Part II) present their 
version of the Elrod algorithm in the context of 
gas dynamics. Borrowing concepts used in model- 
ing transonic flow through a duct, an artificial 
viscosity term is introduced that implements the 
"type differencing" employed. The two schemes 
are equivalent, but the artificial viscosity 
plays an opposite role to Elrod's switch func- 
tion. Their formulation was instructive and not 
unlike the modeling viewpoint of Pan and Frene 
(see discussion on Ranki ne-Hugoniot conditions 
in Part II), although Pan and Frene's treatment 
is more general in that it accomodates any one 
of three types of fluid transport within the 
cavitation boundary. Frene, Launay, and Pan 
illustrate the use of the equivalent Rankine- 
Hugoniot conditions for a centrally fed journal 
bearing using short bearing theory. The analogy 
drawn between the shock wave and the cavitation 
interface is constructive since it forms a basis 
of communication and progress in either field 
can be shared. 

As our models become more complex, the 
search for simpler more efficient programing 
techniques becomes essential. Tichy and Ku 
(Part II) describe a solution procedure for a 
squeeze-film damper bearing with two circumfer- 
ential oil supply grooves at the bearing ends 
undergoing ci rcul ar-centered whirl. In this 
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procedure a large amount of numerical wor< is 
saved and not too much accuracy is lost. Numer- 
ical results for centerline pressure and cavity 
shape using finite-element methods compared 
favorably. 

This brings us to our last topic of discus- 
sion, employing finite-element methods in lubri- 
cation problems in which cavitation occurs. 
Because of the ease in programing and adaptabil- 
ity to the high-speed super computers, finite- 
difference schemes are most often used. How- 
ever, from the standpoint of handling irregular 
boundaries, the finite-element method is well 
suited to dealing with cavitation. In Part II 
Booker discusses several mathematical models 
commonly used to incorporate cavitation in the 
finite-element analysis of lubrication. The 
models are classic in that they are not flow- 
preserving. Until recently, this has been an 
obstacle to achieving a physically realistic 
solution to transient problems. Bayada (Part 
II) discusses a two-dimensional analysis using 
FEM in which mass conservation is treated in a 
way that is consistent with that proposed by 
Elrod (14) or the equivalent Ranki ne-Hugoniot 
condition at the interface. The procedure was 
applied to a seal (46) using rectangular finite 


elements. The discrepancy between input mass 
flow and output mass flow never exceeded 
4 percent in each period. As of this writing, 
Booker (59) has achieved a flow-preserving FEM 
scheme applied to journal bearings in three 
dimensions . 

To sum up, the experimenter must continue 
to help perceive the order in a very complex 
phenomena. This can be done by extensive clas- 
sification and clarification of unresolved 
issues using carefully planned experiments. The 
numerical analyst find themselves on a tortuous 
journey to realism in which all of the effects 
of temperature, misalignment, deflections, local 
deformations, surface energy, and lubricant rhe- 
ology, to name only a few, must eventually be 
included in the computational models. To make 
this possible, the analysts must bring to bear 
the latest numerical advances, keeping abreast 
of new developments in related fields. Algo- 
rithms that incorporate all of the latest tech- 
niques in finite differencing as well as finite 
elements to shorten computational times are 
needed. In addition, these same algorithms must 
lend themselves to employing the latest computer 
architecture, i.e., pipeline and parallel 
processing. 
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